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DISCUSSION 

This  iinal  report  covers  the  period  of  July  1979 
through  April  1980.  and  describes  the  work  performed  under 
contract  N00173-79-C-0450 .  titled  "Radiation  Sources  at 
Electron  Cyclotron  Harmonic  Frequencies." 

This  work  covers  an  initial  investigation  of  gyro- 
tron  devices  operating  in  cylindrical  geometry.  Specific 
topics  include  an  analysis  of  oscillations  in  a  gyrotron 
travelling  wave  amplifier,  the  study  of  the  effects  of 
velocity  spread  and  wall  resistivity  on  gain  and  bandwidth 
in  a  gyrotron  travelling  wave  amplifier,  an  investigation 
of  mode  competition  in  a  cvlindrical  gvrotron  cavity,  and 

V 

finally,  a  numerical  cold  testing  of  cavity  modes  in 
cvlindrical  structures.  _ . 

1.1  Analysis  of  Oscillations  in  the  Gyrotron  Travelling 
Wave  Amplifier 

This  work  was  fully  reported  in  NRL  Memo  Report 
z  4303.  which  is  included  here  as  Appendix  A  to  this  report. 

It  was  found  that  the  convective  instability  which  forms  the 
basis  of  the  gyrotron  travelling  wave  amplifier  may  become 
absolute  (non-convecti ve )  at  a  sufficiently  high  current  level. 
This  results  in  oscillation  rather  than  amplification.  This 
threshold  current  for  the  transition  depends  sensitively  on 
the  applied  magnetic  field.  The  axial  wavelength  and  the 
characteristic  frequency  of  oscillation  at  the  onset  of 
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^  t , 


absolute  instability  are  given.  It  is  round  that  m  »s.er.i  u:: 
spread  has  virtually  no  effect  on  the  threshold  current. 

A  small  amount  of  resistive  wall  loss,  however .  raises  tin 
threshold  current  significantly .  Oscillations  due  to 
part  ieal  reflections  at  the  ends  of  the  system,  are  als> 
examined.  Preliminary  experimental  results  on  both  type 
of  oscillations  are  reported  and  are  found  to  be  in  agree¬ 
ment  with  the  theory. 

1.2  Effects  of  Velocity  Spread  and  Wall  Resistivity  on 
the  Gain  and  Bandwidth  of  the  Gyrotron  Travelling 
Wave  Amplifier 

This  work  was  fully  reported  in  XRL  Memo  Report 
4304.  which  is  included  as  Appendix  B.  to  this  report. 

In  this  report  the  small  signal  gain  bandwidth  of  the 
gyrotron  travelling  wave  amplifier  tgyro-TWA)  operating 
at  the  fundamental  cyclotron  harmonic  is  examined. 

The  analytical  and  numerical  studies  focus  on  the 
effects  of  velocity  spread  in  the  electron  beam  and  of  a 
distributed  wall  resistivity  in  the  waveguide.  It  is 
found  that  wall  resistivity  reduces  the  forward  gain  of 
the  amplifier  only  by  an  amount  approximately  equal  to 
1/3  of  the  corresponding  cold  tube  loss.  Significant 
increase  in  bandwidth  may  result,  under  certain  conditions, 
from  the  presence  of  wall  resistivity.  A  moderate  amount 
(5-10^)  of  velocity  spread  does  not  reduce  the  peak  gain 
significantly  in  general,  but  may  reduce  the  bandwidth  by 
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an  amount  depending  on  t  he  applied  magnetic  iit-ai.  based 
on  c- nisi  derat  i<>ns  oi  such  factors  as  power,  stab:  -  itv . 

^ra  i  n .  bandwidth,  and  eificiency.  we  tentatively  conclude 
that  the  optimum  operating  parameters  for  the  gyro-TWA 
would  be  a  beam  current  on  the  order  of  10  amps,  a  wall 
resistivity  with  a  skin  depth  or.  the  order  oi  5  x  10 
of  the  wall  radius,  and  a  magnetic  field  about  2  percent 
below  the  grazing  value. 
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».o  i/omtH-t  i  p.s:  Moat-  Analysis 

A  study  has  been  made  of  mode  competition  in  a  cylindrical 
gyrotron  cavity.  The  examined  model  is  as  follows:  5efore 
entering  the  cavity,  the  electron  beam  consists  of  equally 
spaced  rings  with  radius  r^,  equal  to  the  Larmor  radius. 

Their  guiding  centers  lie  at  a  distance  ro  from  the  axis  of 
the  cylindrical  cavity.  Each  ring  carries  a  discrete  number 
of  electrons  (usually  16  to  32  electrons).  Initially,  all 
the  electrons  have  the  same  axial  and  normal  velocity.  The 
mode  frequencies  and  spatial  structure  are  unperturbed  by 
the  presence  of  the  beam,  since  the  beam  is  tenuous.  Space 
charge  effects  are  neglected.  Therefore,  the  equations  of 
motion  of  each  electron  are  independent  of  the  other  electrons. 

The  efficiency  for  each  mode  is  computed  by  solving  the 
equations  of  motion  for  all  the  electrons  on  the  rings.  In 
the  nonlinear  regime,  each  efficiency  depends  on  all  the 
amplitudes  of  the  modes  present  in  the  cavity.  It  was  shown 
that  the  efficiencies  converge  to  a  fixed  value  as  the  number 
of  rings  increases  and  all  other  parameters  remain  constant. 

In  each  case  studied,  a  sufficient  number  of  rings  was  chosen 
to  ascertain  convergence. 

For  given  beam  power,  the  amplitudes  of  the  modes  present 
in  the  cavity  were  computed  by  solving  a  nonlinear  set  of 


algebraic  equations  which  express  the  balance  of  power  for 
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each  mode  when  the  system  is  at  a  steady  state.  In  this 
case,  the  power  deposited  to  each  mode  by  the  electron  beam 
is  equal  to  the  power  radiated  by  the  mode  under  consideration. 
The  following  basic  assumption  for  mode  excitation  was  made: 

For  fixed  beam  voltage,  as  the  beam  current  (or  beam  power) 
is  increased  from  zero  to  some  value,  the  system  traces  a 
sequence  of  steady  states  (as  defined  above) .  The  mode  with 
the  lowest  threshold  is  excited  first.  As  the  beam  current 
is  increased  beyond  the  lowest  threshold  value,  a  new  mode 
is  excited  if  this  is  allowed  by  the  presence  of  the*  other 
modes  in  the  cavity  (i.e.  if  the  set  of  nonlinear  algebraic 
equations,  mentioned  above,  does  have  a  solution).  This  is 
a  valid  assumption  since  the  rising  or  decay  time  of  the 
fields  is  much  smaller  than  the  rising  time  of  the  beam  current. 
The  basic  results  are  presented  as  plots  which  express  the 
functional  dependence  of  the  efficiency  and  radiated  power 
of  each  mode  versus  beam  power. 

For  given  beam  power,  a  nonlinear  algebraic  system  may 
have  more  than  a  single  real  solution  and  only  one  of  these 
solutions  is  physically  acceptable.  A  different  approach  to 
the  one  mentioned  above  (namely,  the  solution  of  the  nonlinear 
algebraic  system)  has  been  devised  which  provides  the  necessary 
information  in  selecting  the  physically  acceptable  solution 
for  given  beam  power.  The  information  is  contained  in  contour 


plots  of  certain  functions  which  depend  on  the  amplitudes 
of  the  modes.  These  functions  depend  also  on  the  parameters 
of  the  electron  beam  and  the  cavity  but  they  are  independent 
of  the  Q's  of  the  modes.  Therefore,  the  same  contour  plots 
can  be  used  for  another  cavity  with  different  Q's  for  each 
mode.  This  is  an  advantage  of  this  approach  over  that  of 
solving  the  nonlinear  set  of  algebraic  equations  which  depend 
on  the  specific  Q  values  of  the  modes  for  the  cavity  under 
consideration.  The  contour  plots  have  been  a  valuable  tool 
in  the  determination  of  the  physically  acceptable  solution 
for  given  beam  power,  particularly  in  cases  where  a  jump  of 
the  radiated  power  occurs  for  a  particular  mode  as  the  beam 
power  slightly  increases. 

Two  cases  were  examined : 

i)  TEqh'  teo12  modes  Present  in  the  cavity  and  the 
first  cyclotron  harmonic  is  excited, 
ii)  TE021'  TE022  modes  present  in  the  cavity  and  the 
second  cyclotron  harmonic  is  excited. 

The  preliminary  conclusions  are  as  follows: 

a)  The  Q's  for  each  mode  are  important  in  the 
determination  of  mode  competition. 

b)  For  modes  with  azimuthal  symmetry,  there  is 
indication  that,  as  the  beam  power  builds  up,  the 
mode  with  axial  eigennumber  1=1  will  eventually 
exclude  all  other  higher  order  modes  in  the  cavity. 
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There  is  indication  that  a  node  with  infinite 
threshold  bean  power  may  be  excited  through  the 
excitation  of  other  modes  first  and  subsequent 
mode  conversion. 

and  results  are  presented  in  Appendix  C . 
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and  present  some  sample  results. 

Experimental  void  testing  oi  cavity  structures  t  c 
line  both  ei  geni  requenc  i  es  and  mode  .patterns  is  a  diffi¬ 
cult  problem.  A  numerical  solution  i>  tree  oi  practical 
hardware  problems .  and  can  indeed  give  a  better  answer, 
since  the  cavity  is  in  no  way  perturbed. 

The  section  oi  MASK  relevant  to  this  problem  i s 
the  Maxwell  field  propagator.  If  a  point  disturbance,  a 
delta  function  in  both  time  and  space,  is  introduced  intc 
the  cavity,  the  code  will  properly  evolve  the  resulting 
fields  in  time.  Because  of  the  completeness  Theorem, 
the  initial  delta  function  disturbance  must  contain  all 
the  eigenmodes  of  the  structure. 

As  the  time  proceeds,  the  code  saves  a  complete 
time  history  of  both  the  electric  and  magnetic  fields, 
produced  by  the  disturbance,  at  a  number  of  spatial 
locations.  This  information  can  be  Fourier  analyzed  to 
find  what  frequencies  exist  at  each  location.  This  are 
the  eigenmodes  of  the  cavity. 
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Finally,  t  he  st  ructur*  <.  an  t  •«-  Uri  v*-r.  a:  •a/:. 

:  reiiut-ii ry  :  In:  erest  ler  a  nun.:.-:'  .  :  period.-  .  Th-  r~su  *.  - 
inn  :  it- id  patterns  nan  be  used  to  ider.t  i  !  y  t  h<-  made. 

As  an  example,  a  simple  structure,  cam.- i  s:  inn  <•: 
a  rich:  circular  cvlinuer  was  driver,  by  imposing  a  noise 
puis*.-  at  one  spatial  location,  i  or  m*-  t  i  me  step.  The 
••ode  was  then  allowed  to  run  for  a  number  el  cavity 
transit  times,  with  field  information  saved  at  several 
spatial  locations.  The  use  of  more  than  one  spatial 
location  ensures  that  the  accidental  placement  of  one 
sample  point  at  a  mode  null,  will  not  cause  the  cold 
test  to  miss  a  particular  frequency. 

When  a  Fourier  transform  is  performed  lor  the  data 
stored  at  each  spatial  point,  a  cavity  spectrum  results. 

A  sample  of  such  a  spectrum  is  shown  in  Figure  1.  As 
can  be  seen  from  the  figure,  a  number  of  modes  exist  for 
this  cavity.  Each  peak  can  be  identified  by  driving  the 
cavity  at  that  frequency. 

A  number  of  point  "antennas"  are  placed  into  the 
cavity,  and  they  are  driven  at  the  appropriate  frequency  for 
a  number  of  cycles.  The  TEq^  mode  was  driven,  and  the 
resulting  field  plots  are  shown  in  Figure  2.  In  each 
plot,  the  figure  must  be  rotated  about  the  center  line 
(as  marked)  to  obtain  a  three  dimensional  figure.  Each 
plot  is  a  contour  plot  of  a  particular  field  component. 
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In  summary. this  technique  is  a  straight  forward 
tnethod  tor  determining:  the  resonant  frequency  and  mode 
patterns  lor  arbitrary  cavities.  For  a  simple  cavity, 
as  in  the  example,  field  patterns  identical  to  text  book 
figures  result.  The  real  power  of  the  technique  is  that 
any  arhitrarv  con i i pur a t i on  can  be  tested  with  equal  east 
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Figure  1. 

Cavity  spectrum  of  right  cylindrical 
cavity . 
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ANALYSIS  OF  OSCILLATIONS  IN  THE  GYROTRON 
TRAVELLING  WAVE  AMPLIFIER 


I.  INTRODUCTION 

Current  experiments  performed  at  NRL  on  the  gvrotron  travelling  wave  amplifier  (gyro-TWA) 
have  yielded  record  power  levels  at  millimeter  wavelengths.1 1:1  However,  the  presence  of  unwanted 
oscillations  has  thus  far  prevented  the  operations  of  the  amplifier  in  the  optimum  regime'2'.  Significant 
improvement  in  performance  is  anticipated  with  the  suppression  of  these  unwanted  oscillations.  The 
main  purpose  of  the  present  work  is  to  investigate  the  possible  causes  of  these  oscillations  and  thereby 
to  suggest  ways  for  their  suppression. 

Oscillations  are  generally  caused  by  either  an  external  or  an  internal  feedback  process.  The  exter¬ 
nal  feedback  is  provided  by  reflections  at  both  ends  of  the  interaction  region.  If  the  amplitude  of  the 
reflected  signal  reaches  a  certain  level  that  the  loop  gain  exceeds  unity,  the  wave  amplification  processes 
become  regenerative  and  oscillation  consequently  takes  place.  The  internal  feedback  is  a  result  of  the 
dispersiveness  of  the  unstable  medium.  Under  certain  conditions,  the  wave  may  grow  locally  without 
propagating  axially  out  of  the  system.  As  a  result,  large  amplitude  waves  can  simply  grow  from  noise 
level  perturbations.  Unlike  those  caused  by  end  reflections,  oscillations  produced  by  internal  feedback 
processes  received  comparatively  little  attention.  Their  considerations  have  special  importance  for  the 
gyro-TWA  in  that  the  optimum  operating  regime  of  a  gyro-TWA  happens  to  border  on  the  transition 
from  the  convective  instability  (amplification)  to  the  absolute  instability  (oscillation).  Careful  study 
of  the  internal  feedback  processes  can  therefore  provide  a  guide  for  acheiving  optimum  operation  of  the 
gyro-TWA  without  the  excitation  of  absolute  instabilities. 

Experimentally,  there  is  convincing  evidence  that  both  types  of  oscillations  are  excited  in  the 
gyro-TWA  operated  at  NRL.  The  two  types  of  oscillations  have  different  characteristics.  Toward  the 
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end  of  this  paper,  we  shall  briefly  describe  the  experimental  results  in  the  light  of  the  present  theoreti¬ 
cal  study.  The  details  of  these  experiments  have  been  reported  elsewhere. 

We  shall  first  (Sec.  II]  present  a  quantitative  study  of  the  threshold  current  above  which  an  abso¬ 
lute  instability  would  be  excited  as  a  result  of  internal  feedback  processes.  The  influence  of  velocity 
spread  and  resistive  wall  loss  on  the  threshold  current  will  also  be  examined.  We  next  consider  oscilla¬ 
tions  caused  by  external  reflections  in  Section  III.  where  experimental  observations  of  both  types  of 
oscillations  are  reported. 


II.  INTERNAL  FEEDBACK  MECHANISMS-CONDITION  FOR 
THE  ONSET  OF  ABSOLUTE  INSTABILITY 

When  operated  near  the  "grazing  condition,”  the  mildly  relativistic  electron  beam  of  the  gvro- 
TWA  interacts  most  effectively  with  the  waveguide. 14  51  This  "grazing  condition"  is  achieved  when  the 
longitudinal  velocity  of  the  beam  is  equal  to  the  group  velocity  of  the  forward  wave  of  the  waveguide 
mode.  If  the  beam-wave  coupling  (determined  by  the  current  strength  and  magnetic  field  match,  etc.) 
is  sufficiently  weak,  the  beam  will  only  couple  to  the  forward  waveguide  mode.  The  gyro-TWA  would 
then  operate  according  to  what  is  expected  of  a  convective  instability. iV  However,  if  the  coupling  is 
sufficiently  strong,  the  backward  wave  of  the  guide  mode  would  be  excited.  The  excitation  of  these 
backward  waves  might  then  provide  an  internal  feedback  mechanism,  resulting  in  an  absolute  instability 
which  would  manifest  itself  as  a  natural  oscillation  of  the  system.  (Under  these  conditions,  care  must 
be  exercised  to  intepret  the  gain  curve  of  the  "amplifier".'4') 

The  critical  current  which  marks  the  transition  to  absolute  instability  depends  sensitively  on  the 
magnetic  field.  It  is  obtained  by  using  the  stability  criterion  of  Briggs  and  Bers. ,JI  The  analysis  is  based 
on  the  assumption  that  the  system  is  infinitely  long.  The  characteristic  frequency  and  wavenumber  at 
the  transition  based  on  such  an  analysis  would,  however,  shed  light  on  the  required  length  (starting 
length)  of  the  system  if  it  is  to  avoid  these  natural  oscillations. 
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The  result  ts  summarized  as  follows.  If  both  velocity  spread  and  wall  resistivity  are  absent,  and  if 
the  applied  magnetic  field  is  maintained  near  the  grazing  condition,  a  convective  instability  sets  in  at  a 
relatively  low  current.  This  convective  instability  becomes  an  absolute  one  if  the  current  (beam 
strength)  exceeds  a  certain  critical  value.  This  critical  value  is  given  analytically  [cf.  Eq.  (10)]  in  terms 
of  the  waveguide  dimension,  electron  parallel  and  perpendicular  velocities,  and  the  applied  magnetic 
field.  At  the  onset  of  the  absolute  instability,  the  oscillation  frequency  and  the  axial  wavenumber  are 
also  given  analytically.  It  is  found  that  typically  the  oscillation  frequencies  are  almost  identical  (but 
slightly  less  than)  with  the  cut-off  frequency  of  the  waveguide.  We  found  that  velocity  spread  in  the 
beam  has  virtually  no  effect  on  the  threshold  current  for  the  transition  to  absolute  instability.  The 
presence  of  a  small  amount  of  distributed  wall  resistivity  would,  however,  raise  the  threshold  current 
significantly.  It  can  be  shown  that  such  a  distributed  wall  loss  does  not  significantly  reduce  the  gain  of 
the  amplifier.1”  In  this  regard,  the  presence  of  wall  resistivity  offers  the  best  hope  to  render  the 
amplifier  stable  at  high  current  level. 

In  Section  Il-a,  we  present  the  details  of  the  analysis  for  the  case  of  a  cold  beam  situated  in  a  loss¬ 
less  waveguide.  The  method  developed  there  has  been  extended  to  the  studies  of  a  beam  with  finite 

temperature  in  Section  Il-b.  Finally,  in  Section  II-c,  we  present  the  results  for  a  lossy  waveguide. 

* 

a.  Cold  Beam,  Lossless  Waveguide: 

The  typical  configuration  of  a  gyro-TWA  (Fig.  1)  consists  of  an  annular  electron  beam  propagat¬ 
ing  inside  a  waveguide  of  circular  cross  section  of  radius  rw.  The  electrons,  guided  by  a  uniform  mag¬ 
netic  field  ( B0  e;),  move  along  helical  trajectories.  Ideally,  all  electrons  have  the  same  perpendicular 
velocity  v10,  and  the  same  parallel  velocity  (v;0),  with  their  guiding  centers  uniformly  distributed  on  a 
surface  of  constant  radius  r0  (Fig.  1).  The  dispersion  relation  for  the  interaction  between  this  cold 
beam  and  the  TEm„  waveguide  mode  (m,  n  being,  respectively,  the  azimuth'.!  and  radial  mode 
numbers)  at  the  s-th  cyclotron  harmonic  frequency  has  been  derived  in  Ref.  (5).  It  reads 
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,  ,  -  ,  .  -  -  4i >c:  ~  k.:c:) 

<u‘  -  k.'C  -  k+rC‘  “  - ; -  - : - 

>or;Kmr.  lu-  -  k.\;0-  sH.i- 

(o;  -  k.y;o)  Qv T  ( kmrr,j,  k„ri0 ) 
to-*.V.0-jfl. 

where  to  is  the  frequency  of  the  mode.  k.  is  the  axial  wavenumber.  y0  -  tl  -  \_:0  c:  -  \:l/  r:>~ 

(3.o  ”  vi(/c.  n,  -  eflo/To  we.  rip  -  v  ^n..  kmn  -  X„Jrm%  Xmn  is  the  n-th  root  of  J‘ mKx'  -  0  J„ (*•  is 
the  Bessel  function  of  order  m.  J’m(x)  -  dJm(x)/dx,  v  =  Se2.1  me2  is  a  dimensionless  beam  density 
parameter,  .V  is  the  number  of  electrons  per  unit  axial  length,  and  the  functions  Km„.  Him.  Qim  are 
defined  as  follows: 


Kmn  -  J2{Xmn)  U  - 


HsmUy)  -  y,_mU)y's(.v) 


Qsm(x.y)  -  2/fjm(x>')  +  y  J-.  m  (x)J\{y)J"s(y)  +  —  J2.  m  .  ,(x) 


■  AO-)/’,. ,0-)  -  y  U)A0')A*i  O')]. 

Instability  occurs  at  wave  frequency  and  axial  wavenumber  such  that  w:  -  k?c 2  -  k2„c:  *sOand 
to  -  k. v;0  -  sfl.  =  0.  Hence  the  second  term  on  the  right  hand  member  of  Eq.  (1)  can  usually  be 
ignored  in  comparison  with  the  first  term.  We  can  further  approximate  tu:  -  k.2c2  by  k2„c2  in  the  right 
hand  side  of  (1).  Thus,  Eq.  (1)  is  simplified  to 

(to:  -  k?c2  -  k2„c2)  (to  -  k.v:0—  jfl  r)2  - - *Tr~  P''o  Am-  <2) 

y  Xmn 

Henceforth,  we  shall  concentrate  on  the  stability  analysis  according  to  Eq.  (2). 


It  is  convenient  to  normalize  the  frequency  w  with  respect  to  the  cutoff  frequency  a»c  s  A.„„cand 
k.  to  km„.  Then  Eq.  (2)  is  written  in  dimensionless  form  as 

(S2-  P-  1)(S  -  ik/3,,-  b):-  -e,  (3) 

where  S  s  w/tof;  k  s  k./k„„;  /9«  =  v.0 /c;  and  the  dimensionless  parameters  band  c  are  defined  by 


b  s  s  Cl  c/a>c 
yo'X.o't 
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Note  that  b  is  proportional  to  the  magnetic  field  and  t  is  a  measure  of  the  beam  strength.  Typically,  « 


is  a  very  small  number,  on  tne  order  of  1(T5  tor  even  less'.  As  an  example12',  take  m  »■  0.  s  «  n  - 
l.  “  -4.  &  “  .266,  <u,/ 2 ~  “  34.3  GHz.  r .  —  .533  cm.  Then 

<  -  3.78  x  l(T6  /  14c) 

where  I  is  the  beam  current  in  units  of  ampere. 

In  the  limit  of  zero  beam  strength  U  -  0)  the  waveguide  modes 

-  P  -  1  -  0,  (5a) 

and  the  beam  mode 

to  -  k(3r,  -  b  -  0.  (5b) 


are  decoupled  [cf.  Fig.  (2)].  The  coupling  of  these  two  modes  is  strongest  when  the  dispersion  curve  of 
the  beam  mode  is  tangent  to  the  waveguide  mode  [cf.  Fig.  (2)).  This  occurs  when  the  magnetic  field  is 
adjusted  so  that  b  -  b0,  where 

b0={  l-/3,;)1/:.  (6) 

In  the  following,  we  shall  assume  that  /3,  is  fixed,  and  that  b  is  in  the  vicinity  of  b0.  We  shall  then 

study  the  properties  of  the  dispersion  relation  (3)  as  «.  the  beam  strength,  is  varied.  This  way  of  fixing 

parameters  is  suitable  for  comparison  with  experiments." 

We  state  here  the  principal  results.  The  details  of  derivations  of  these  results  are  given  in  the 
Appendices. 

(i)  If  b  ^  b0 ,  there  is  an  instability  whenever  «  >  0.  This  instability  can  either  be  convective  or 
absolute.  These  statements  are  proven  in  Appendix  A. 

(ii)  If  b  -  b0  (i.e..  exact  grazing  condition),  there  is  a  convective  instability  if  0  <  e  <  e cs.  This  ins¬ 
tability  becomes  an  absolute  one  when  t  >  t  C!  where  e  ct  is  defined  by 


(27) w4  y/jr,  kS(, 

(7) 

ky  “  (toJlt  —  ip)  /  4  /3ii, 

(8) 

(60  +  6v2  0 
*  1  +  8  fit 

(9) 
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Furthermore.  at  the  onset  of  the  absolute  instability,  (i.e..  «  “  t,f),  the  characteristic  wavenumber  k:i 
and  frequency  Zsg  are  given  by  Eqs.  (8'  and  (9)  respectively.  Note  that  the  wavenumber  at  the  onset 
of  absolute  instability  is  purely  real.  This  general  result  is  demonstrated  in  Appendix  B,  The  proof  of 
the  rest  of  the  statements  in  this  paragraph  is  given  in  Appendix  C. 

(iii)  For  b  —  b0.  the  onset  of  the  absolute  instability  occurs  when  <  —  where  t.  is  defined  by 


(27) 

**  “  40~  (Sj  ~  b) 1 


(10) 

(ID 


b  -r  [8/3  2(1  —  b~)  -r  64fl.  ]  * 

to, - e -  (12) 

1  -  8/3- 

Again,  ki  and  5,  are  respectively  the  wavenumber  and  frequency  of  oscillation  at  the  onset  of  the  abso¬ 
lute  instability.  It  is  clear  that  Eqs.  (10M12)  reduce  to  (7)-(9)  when  b  -  b0.  i.e..  when  the  gyro-TW.A 
is  operated  exactly  at  the  grazing  condition. 


The  numerical  values  of  the  critical  beam  strength,  as  well  as  the  corresponding  frequency  of 
oscillations  and  wavenumber,  can  readily  be  computed  for  various  values  of  0V  and  b  from  Eqs.  (7>- 
(12).  Listed  in  Table  1  are  the  values  generated  from  Eqs.  (7)- (9)  which  physically  correspond  to  the 
exact  grazing  condition.  Note  from  the  last  column  of  this  table  that  the  threshold  beam  current 
(measured  by  tsf)  increases  with  )3„.  This  result  can  readily  be  threshold  beam  current  (measured  by 
<w)  increases  with  /3„.  This  result  can  readily  be  understood  by  noting  that  at  higher  axial  beam  velo¬ 
city.  disturbances  conveci  away  at  a  higher  rate.  To  couple  the  forward  propagation  of  disturbances 
with  the  backward  waveguide  mode,  a  condition  for  the  excitation  of  absolute  instabilities'-1,  a  stronger 
current  is  required.  In  Table  2.  we  fix  /3.  at  .266.  We  then  calculate  the  threshold  current  at  various 
values  of  magnetic  field  strength.  We  note  that  if  the  magnetic  field  is  below  the  grazing  value,  a 
higher  beam  current  is  necessary  in  order  to  excite  the  absolute  instability,  because  the  beam  mode  and 
the  waveguide  mode  are  away  from  perfect  synchronism.  If  the  magnetic  field  is  above  the  grazing 
value,  the  unstable  spectrum  readily  extends  to  negative  values  of  k.  and  the  absolute  instability  is 
easily  excited.  In  the  event  that  the  beam  mode  (Fig.  2)  intersects  with  the  backward  waveguide  mode 
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'.he  absolute  lnstaoility  sets  ir.  as  long  as  the  beam  current  is  non  zero.  In  the  last  two  columns  of 
Table  2.  we  have  assumed  that  the  threshold  current  /  is  related  to  «  by  Eq.  <4c>  Note  further  that  tne 
oscillation  frequency  of  these  oscillations  is  slightly  less  than  the  cutoff  frequency  of  the  waveguide. 

In  the  current  gyro-TWA  experiments  at  NRL,i:'.  £  “  .266.  cu./2 it  -  34.3  GHz  and  kmr  - 
2.18  :m~'  for  the  TE01  mode.  If  the  magnetic  field  is  maintained  at  the  grazing  condition  (6/60  “  1) 
then  an  oscillation  at  a  frequency  of  34.26  GHz  is  predicted  to  occur  when  the  beam  current  is  raised 
beyond  .588  amps  as  a  result  of  absolute  instability  [cf..  Row  (S'  in  Table  II].  Experimentally,  an  oscil¬ 
lations  at  34.27  GHz  was  indeed  observed'1'  at  a  beam  current  of  1  amp  as  the  magnetic  field  is  raised 
to  the  vicinity  of  the  grazing  condition.  The  actual  threshold  current  observed  should  be  somewhat 
higher  than  the  predicted  value  as  the  latter  was  made  for  an  infinitely  long  waveguide  whereas  in  real¬ 
ity  a  finite  length  waveguide  was  used  and  coupling  out  the  ends  was  not  included  in  the  present  calcu¬ 
lation.  We  shall  discuss  in  greater  detail  the  experimental  results  in  Section  III. 

b.  Warm  Beam.  Lossless  Waveguide: 

We  have  studied  the  effects  of  the  velocity  spread  and  energy  spread  on  the  threshold  current 
above  which  an  absolute  instability  is  excited.  In  the  NRL  gyro-TWA  experiments,  the  energy  spread 
is  minimal,  of  order  1%  or  less.  The  spread  in  the  longitudinal  velocity  is  considerably  higher,  on  the 
order  of  10  per  cent.  We  found  that,  even  if  the  velocity  spread  is  as  high  as  20  percent,  it  has  a  negli¬ 
gible  influence  on  the  threshold  current.  In  other  words,  momentum  spread  cannot  stabilize  the  abso¬ 
lute  instability.  This  seemingly  surprising  result  can  be  understood  as  follows.  The  presence  of  an 
absolute  instability  requires  the  backward  waveguide  mode  to  be  coupled  to  the  forward  beam  mode13'. 
The  axial  wavenumber  k.  is  thus  constraint  to  a  value  in  the  vicinity  of  zero  at  the  transition.  [See  also 
column  3  of  Table  II].  Now,  the  velocity  spread  <dv:>  modifies  the  dispersion  relation  by  an  amount 
on  the  order  of  ( k .  <Av.>  Vwf:,  which  is  an  exceedingly  small  quantity  at  the  transition  to  the  absolute 
instability  by  virtue  of  small  value  of  k.,  leaving  the  threshold  current  calculated  in  the  previous  sub¬ 
section  essentially  unchanged. 
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The  above  conclusions  stem  from  our  analysis  of  the  dispersion  relation  which  includes  the  effect 
of  velocity  spread  and  energy  spread,  measured  respectively  by  <Av;>  v..  and  <Ay  >  -y,,.  Assuming 
that  these  quantities  are  small,  we  obtain  the  dispersion  relationship 


* 

A 

_ 

_  <  Av.> 

_  <Av> 

* 

k'  -  I ) 

—  a) 3  —  b  )*  — 

kB  - ; — 

— 

<JJ  '  — 

V;0 

>c 

This  equation  can  readily  be  compared  with  Eq.  (3).  where  the  symbols  are  defined.  Equation  (13)  can 
be  deduced  from  the  one  obtained  by  Unm19'  who  applied  the  idealized  ’water-bag'  electron  beam  dis¬ 
tribution  function  to  the  formulation  developed  earlier.'10’ 

Applying  the  techniques  developed  in  the  Appendices,  we  found  that  a  velocity  spread  less  than 
twenty  percent  has  virtually  no  effect  on  the  transition  current,  for  reasons  explained  above.  An  energy 
spread  less  than  half  a  percent  actually  slightly  lowers  the  threshold  current. 

It  should  be  emphasized  that,  while  the  momentum  spread  is  not  effective  in  suppressing  the 
absolute  instability,  it  can  reduce  the  gain  and  bandwidth  of  the  amplifier  significantly'11.  Once  the  fre¬ 
quency  of  the  input  signal  is  above  the  cutoff  frequency  of  the  waveguide,  k.  will  increase  significantly 
(in  comparison  with  a  value  close  to  zero).  The  velocity  spread  will  then  have  a  significant  effect  on 
the  gain  [cf.  Eq.  (13)].  This  observation  is  corroborated  by  the  numerical  studies  on  the  gain  and 
bandwidth  of  the  amplifier  reported  in  Refs.  (1)  and  (7). 

c.  Cold  Beam.  Lossy  Waveguide: 

An  attempt  to  stabilize  the  gyro-TWA  is  to  coat  the  inner  wall  of  the  metallic  waveguide  with  a 
layer  of  material  of  finite  electrical  conductivity,  such  as  graphite.'8’  A  detail  study  of  the  effect  of  a 
resistive  wall  on  the  gain  of  the  gyro-TWA  is  beyond  the  scope  of  this  paper.  Here  we  concentrate 
only  on  the  dependence  of  threshold  current  on  wall  resistivity.  It  is  found  that  even  a  small  amount 
of  wall  loss  increases  the  threshold  current  substantially. 

Let  8  be  the  skin  depth  of  the  wall  material  of  the  waveguide.  Typically.  8/r.  «  1  where  r,  is 
the  wall  radius.  In  the  limit  of  small  8/r.,  the  dispersion  relation  reads'  ’ 
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which  is  again  to  be  compared  with  Eq.  <3)  In  obtaining  Eq.  <  1  -i ) .  we  have  ignored  the  effect  of  velo¬ 
city  spread  as  it  can  hardly  modify  the  threshold  current.  In  passing,  we  remark  that  in  the  resistive 
wall  experiments  to  be  performed  at  NRL18'.  6  r„  s  2.6  x  10~-\  and  that  values  of  6  r,  as  nigh  as  10': 
are  reasonable  in  practice. 

The  analysis  of  Eq.  (14)  for  the  threshold  current  is  more  involved  and  will  not  be  detailed  nere. 
Shown  in  Fig.  (3)  are  the  results.  To  obtain  the  data  in  this  figure,  we  fix  13  —  .266  The  various 
curves  correspond  to  different  magnetic  fields.  They  illustrate  the  threshold  current  as  a  function  of 
hr,  iwall  resistivity).  It  is  seen  that  a  small  amount  of  wall  resistivity  can  significantly  increase  the 
threshold  current.  That  is.  wall  resistivity  may  render  the  gyro-TW’A  stable  against  the  absolute  insta¬ 
bility.  Shown  in  Table  III  are  the  data  for  b /  b0  —  ,98  at  the  transition.  Again,  we  have  assumed  that 
the  beam  current  I  is  related  to  t  by  Eq.  (4c)  in  obtaining  these  results. 

The  increase  of  threshold  current  with  wall  resistivity  may  be  understood  qualitatively  as  follows. 
As  the  frequency  of  oscillation  at  transition  is  below  the  cutoff  frequency  of  the  waveguide,  the  back¬ 
ward  wave  of  the  waveguide  mode  attenuates  in  the  reverse  direction  of  beam  propagation.  This 
attenuation  rate  increases  as  wall  resistivity  (6/ r,)  increases.  For  an  absolute  instability  to  be  excited, 
ihis  attenuation  rate  must  match  the  amplification  rate  of  the  forward  beam  mode.  This  matching 
requirement  is  just  the  alternative  statement  that  a  saddle  point  (i.e..  merging  of  roots  of  k  at  a  single 
frequency)  must  be  present  for  the  onset  of  absolute  instability.'3’  A  higher  beam  current  is  thus 
needed  to  excite  the  absolute  instability  when  wall  resistivity  is  present.  Note  again  that  the  oscillation 
frequency  is  just  below  the  cutoff  frequency  of  the  waveguide  at  the  transition,  and  that  k,  is  complex 
for  reasons  given  above. 

We  point  out  that  the  internal  feedback  mechanism  has  recently  been  theoretically  analyzed  by 
Etlicher  and  Buzzi."1'  However,  the  main  stabilizing  influence  of  wall  resistivity  has  not  been  included 
in  their  calculations.  Their  work,  reported  in  a  conference  abstract"11,  does  not  contain  sufficient  detail 
for  a  direct  comparison  with  our  work. 
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III.  EXTERNAL  FEEDBACK  MECHANISMS  AND  EXPERIMENTAL  OBSERVATIONS  OF 
OSCILLATIONS. 

Treating  the  length  of  the  amplifier  as  infinite  we  study  in  the  previous  section  the  convective  and 
absolute  nature  of  the  cyclotron  maser  instability.  When  an  absolute  instability  first  sets  in.  the  fre¬ 
quency  of  oscillations  is  predicted  to  be  just  slightly  less  than  the  cutoff  frequency  of  the  waveguide. 
The  axial  wavenumber  is  typically  1/30  the  transverse  wavenumber  of  the  waveguide.  Thus,  a  short 
interaction  length  (as  well  as  a  resistive  wall)  may  reduce  the  danger  of  exciting  absolute  instabilities. 

In  reality  the  gyro-TW.A  system  is  finite  in  length  and  an  additional  mechanism  for  oscillation  is 
provided  by  partial  reflections  of  the  amplifying  w'aves  at  the  ends  of  the  system.  At  high  beam 
currents,  the  e-folding  length  of  spatial  amplification  is  very  short,  and  a  slight  mismatch  may  convert 
an  amplifier  into  an  oscillator.  Based  on  this  argument,  one  can  see  that  oscillations  due  to  end 
reflection  will  inevitably  set  in  at  a  high  current  level,  regardless  of  how  small  a  mismatch  there  is  at 
the  ends  of  the  device.  Oscillations  induced  by  end  reflections,  as  well  as  their  suppressions,  are  well- 
known  in  the  technology  of  travelling  wave  tubes.  For  the  present  case  of  the  gyro-TWA.  this  type  of 
oscillations  will  have  frequencies  higher  than  the  cut-off  frequency  of  the  waveguide  for  only  then 
would  the  gain  of  the  amplifier  be  significant.1 11  These  statements  appear  to  have  been  corroborated  by 
the  NRL  experiments  to  be  decribed  below.  Again,  a  resistive  wall  can  (1)  render  the  amplifier  stable 
against  this  type  of  oscillations,  and  at  the  same  time.  (2)  increase  the  bandwidth  of  the  amplifier, 
without  significantly  sacrificing  its  gain.'71 

Instabilities  have  been  routinely  observed<:'8'  during  experimental  operation  of  gvro-TWA 
operating  in  the  T£0i  circular-electric  mode  at  35.1  GHz  at  NRL.  The  physical  design  of  the  gyro-TWA 
uses  a  0.533  cm  radius,  for  a  cutoff  frequency  of  34.30  GHz.  copper  waveguide  21  cm  long.  The 
parameters  of  the  beam  were  designed  to  be  v,,  -  .27  c  and  v:  .40  c.  Typical  linear  (small  signal) 
gains  obtained  were  24  db  at  a  3  ampere  beam  current  and  32  db  at  9  amps. 

Typical  oscillations  occur  at  several  discrete  frequencies  normally  just  above  cutoff  at  34.30  GHz. 

such  as  34.38  GHz  and  34.45  GHz.  The  power  output  of  these  oscillations  is  typically  on  the  order  of 
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;he  maximum  saturated  output  during  stable  amplification.  These  oscillations  can  be  induced  from  a 
staoi>  operating  amplifier  by  increasing  the  beam  current  and/or  the  magnetic  fieid  hence  increasing  the 
gam  These  oscillations  can  be  explained  as  follows. 

The  gyro-TWA  is  operated  near  the  cutoff  frequency  of  the  interaction  waveguide.  The  mechan¬ 
ism  actually  has  considerable  gain  almost  dow’n  to  cutoff.  ‘  Even  though  the  input  and  output  couplers 
were  carefully  designed  for  low  reflections  throughout  the  useful  bandwidth  of  the  gvro-TWA.  match¬ 
ing  becomes  impossible  as  the  cutoff  frequency  is  approached.  Thus,  reflection  spikes  are  obtained  near 
cutoff.  Cold  test  reflection  measurements  indeed  reveal  reflections  at  the  problem  frequencies,  namely, 
at  34.38  GHz  and  34.45  GHz. 

As  the  gyro-TWA  beam  current  or  magnetic  field  is  increased  the  gain  near  cutoff  is  increased. 
The  reflections  constitute  a  feedback  and  when  the  net  loop  gain  exceeds  loop  losses  then  oscillation 
sets  in.  This  is  evidenced  by  the  oscillations  at  34.38  GHz  and  at  34.45  GHz  reported  above. 

Instability  at  just  beiow  empty  waveguide  cutoff  has  also  been  observed  and  is  believed  to  be  the 
absolute  instability  predicted  in  Sec.  II.  The  conditions  are  as  follows.  At  a  beam  current  of  1  amp  and 
cathode  voltage  of  -70  kV,  the  magnetic  field  was  scanned.  The  gain  measured  at  35.13  GHz  peaked  at 
12.7  dB  with  1230  watts  output  (unsaturated)  at  a  magnetic  field  of  12.81  kG.  As  the  magnetic  field 
was  increased  the  gain  (at  35.13  GHz)  dropped  rapidly  to  —0  db  at  13.11  kG.  Self-oscillation  in  the 
TEo\  mode  at  34.45  GHz  set  in  at  a  magnetic  field  of  13.15  kG  at  a  power  output  level  of  8.71  kW. 
This  power  level  is  on  the  order  of  maximum  observed  saturated  amplifying  power  level.  An  observed 
(stable)  maximum  output  power  at  saturation  at  1  amp  beam  current  was  5.75  kW  with  4.3  db  gain.  As 
the  magnetic  field  was  raised  the  oscillations  at  34.45  GHz  continued.  At  13.37  kG  the  34.45  GHz 
oscillation  ceased  and  another  oscillation  at  34.27  GHz  was  observed.  This  oscillation,  slightly  below 
the  TE§ ,  cutoff,  had  an  output  power  of  200  watts.  At  13.41  kG  the  power  output  in  this  oscillation 
increased  to  660  watts. 
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The  experimental  data  reported  here  is  limited.  At  the  time  the  data  was  taken  1  we  were 
unaware  of  potential  absolute  instabilities  or  their  characteristics.  Close  attention  will  be  made  to  abso¬ 
lute  instabilities  in  future  experiments.  However,  we  tentatively  interpret  the  above  data  as  follows. 

As  the  magnetic  field  is  increased  to  the  peak  gain  point  measured  at  35.13  GHz  (corresponding 
to  a  grazing  condition),  the  gain  near  cutoff.  i.e.,  34.30  GHz  to  34.50  GHz.  is  too  low  to  sustain 
reflection  type  oscillations.  As  the  magnetic  field  is  increased  past  the  peak  gain  point  (at  35.13  GHz) 
the  gain  is  still  increasing  near  cutoff  until  sufficient  loop  gain  is  attained  for  oscillation  at  34.45  GHz. 
As  the  magnetic  field  is  increased  beyond  grazing,  the  threshold  current  for  the  absolute  instability 
rapidly  decreases.  When  the  absolute  instability  becomes  dominate,  then  oscillation  occurs  in  that 
mode.  The  calculations  of  the  threshold  indicate  that  the  absolute  instability  should  occur  even  at  graz¬ 
ing.  i.e..  the  calculated  threshold  current  is  0.588  amp.  However,  the  calculation  was  made  for  an 
infinitely  long  waveguide  when  realistically  a  finite  length  waveguide  was  used  and  coupling  out  the 
ends  was  not  included.  Thus  we  would  expect  the  actual  threshold  to  be  somewhat  higher  than  indi- 
cated  in  Table  II.  The  lower  power  level  output  seems  reasonable  since  the  power  is  essentially  coupled 
out  the  end  rather  than  propagated  out. 

We  then  have  several  points  supporting  our  belief  in  the  observation  of  the  absolute  instability, 
(i)  The  oscillation  was  slightly  below  cutoff  in  TE^  mode  exactly  at  the  predicted  value:  (ii)  the  power 
output  level  was  substantially  lower  than  the  routinely  observed  reflective  oscillations  indicating  a 
different  type  of  instability;  and  (iii)  as  predicted,  the  absolute  instability  set  in  as  the  magnetic  field 
was  raised.  Although  the  reflection  type  instability  seems  to  be  the  most  dominant,  especially  at  higher 
current  levels,  this  absolute  instability  could  be  important  in  gyro-TWA  design  with  regard  to  methods 
of  overall  stabilization. 

An  attractive  method  to  suppress  both  types  of  oscillations  is  to  add  loss  to  the  circuit  to  decrease 
the  loop  gain.  This  also  raises  the  threshold  current  for  the  excitation  of  absolute  instabilities.  Further¬ 
more.  calculations'71  show  that  the  forward  gain  in  a  distributed  loss  waveguide  circuit  is  only- 
decreased  by  approximately  one-third  the  cold  tube  loss  in  dB  per  unit  length.  Thus  stability  would  be 
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obtained  with  only  a  slight  decrease  in  forward  gain.  In  addition,  wall  loss  leads  to  bandwidth  improve¬ 
ment. 
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APPENDIX  A 


Presence  of  an  instability  Whenever  b  >  i>0. 


In  this  Appendix,  we  prove  that  there  is  always  an  instability  (which  can  either  be  absolute  or 
convective)  whenever  the  magnetic  field  is  higher  than  that  at  the  grazing  condition.  It  suffices15'  to 
show  that,  whenever  «  >  0  and  b  ^  b0 ,  there  is  always  a  pair  of  complex  roots  of  for  some  real 
values  of  k.  say.  at  k  -  £„/6.  With  £  -  0„/6,  Eq.  (3)  can  be  written  as 

(jr  -  l)(y  -  a)2  «  -  (Al) 

where 

v  2  ub/(b2  +  (A2) 

a  =  (£,:  +•  b •)'  (A3) 

e,  2  €  64/(*;  +  (A4> 

Note  that 


since 


a  >  1 


(A5) 


0,?  +  b:  -  1  +  (i:  -  b$)  2  1  if  b>  bo- 

As  a  function  of  y ,  /(.v)  2  (y:  -  l)(y  -  o):  is  sketched  in  Fig.  4.  As  easily  seen  from  this  figure,  the 
equation  f(y)  —  admits  at  least  a  pair  of  complex  roots  of  y  whenever  «,  >  0.  From  (A2)  and 
(A4).  we  conclude  that  w  has  at  least  a  pair  of  complex  roots  as  long  as  <  >  0. 
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Onset  of  Absolute  Instability  at  Real  Values  of  L 

In  this  Appendix,  we  prove  that  when  an  absolute  instability  first  sets  in.  the  wavenumber  is 
always  real.  This  result  is  useful  because  it  considerably  restricts  our  search  in  wavenumber  space  in 
the  calculation  of  the  critical  beam  strength  that  marks  the  transition  to  absolute  instability. 

At  the  transition,  the  frequency  Z  -  Zs  is  real.  For  this  value  of  Z,.  there  is  a  double  root13'  of 
k.  Below,  we  prove  that  this  double  root  is  real. 

Since  Z  -  Zs  is  real,  and  Eq.  (3)  is  a  fourth  degree  polynomial  of  k  with  real  coefficients,  all  four 
roots  k i.  £3.  kt  of  Eq.  (3)  must  either  be  real  or  appear  as  complex  conjugate  pairs.  We  shall  first 
assume  that  the  above-mentioned  double  root  has  a  nonvanishing  imaginary  part  and  then  show  that  a 
contradiction  results. 


Denote  the  real  and  imaginary  part  of  the  double  root  it,,  k2  to  be  ksr  and  ks,.  If  ksi  ^  0.  the 
other  two  roots  £5,  it*  must  also  be  a  double  root  since  the  roots  occur  in  conjugate  pairs,  in  this  case. 
Eq.  (3)  may  be  written  as 


(k  -  kx)(k  -  k:)(k  -  it 3) (£  -  kj 
“  (k  -  k„-  i  kv)2(k  -  kv  +  i  k„)2 

-  it4  +  P  ( — 4 A:,,)  +  P( 6 k2  +  2 k2) 

+  k  [—Aksr(k2  +  k2)}  +  ( k2  +  k2)2  -  0. 

On  the  other  hand,  equation  (3)  may  also  be  written  as 


it4  +  P  2 


b  - 


+  if 


+  ir  -2  (Z2  -  1) 

b  —  Z, 

-  -  (Zj  -  1) 

b  —  Z  j 

/3. 

/8.. 

(Bl) 


(B2) 
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Since  (Bl)  and  (B2)  are  the  same  equation,  the  coefficients  must  be  the  same.  Thus. 


6k- -  2k- - 


—  4  ks,  “  2(b  -  a»,)//3  , 
b  — *  u> 


-  4 *„(*  =  +  fcs;)  -  -2  (S;  -  1) 
Dividing  (B3)  by  (B5).  we  obtain 

-  (Ss2  -  1)  -  k;r  +  fc;. 


b  —  at, 


—  (Zjs‘  —  1)  —  5 k,‘, 


( B3> 

i). 

(B4) 

b  —  u>  s 

(B5) 

0 

(B6) 

0  as 

;  k-. 

(B7) 

Clearly,  (B7)  contradicts  (B4)  if  k„  ^  0.  We  conclude  then  that  the  assumption  that  k%  =*  0  is  false. 
Hence  at  the  transition  to  absolute  instability,  the  wavenumber  k,  must  be  real. 
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Detailed  Study  of  Absolute  and  Convective  Instability  at  the  Grazing  Condition. 

We  shall  focus  our  attention  to  the  case  where  b  -  b0.  The  dispersion  relation  (3)  reads 

D(u).k)  s  [k:  -  (w:  —  1)1  Ar  —  W-~-  -  -  —  (Cl) 

&  fir 

where,  in  this  Appendix,  we  have  dropped  the  bar  in  w  and  k  for  convenience.  Since  the  beam 
strength  t  enters  only  in  'he  RHS  of  (Cl),  we  shall  examine  the  detailed  property  of  D(u.k).  In  the 
present  study  of  transition  to  absolute  instability,  we  need  only  to  consider  real  values  of  a>  and  k.  [Cf. 
Fig.  51.  Since  the  onset  of  an  absolute  instability  always  occurs  at  a  multiple  root  of  real  k  (Cf.  Appen¬ 
dix  B).  we  need  only  to  focus  our  attention  at  those  values  of  k  in  Figs.  5  for  which  D(w.k)  is  station¬ 
ary.  Moreover,  since  «  ^  0,  only  non-negative  values  of  D(w.k)  are  considered. 

We  first  determine  the  trajectories  of  the  roots  in  the  complex  k  plane  as  the  imaginary  part  of  w 
changes  from  0  to  [Cf.  Fig.  6a.]*  This  consideration  is  necessary  to  distinguish  an  absolute  from 
convective  instability. 131  As  <u  —  —  /  °°,  among  the  four  roots  of  k  in  Eq.  (Cl),  three  roots  would  tend 
to  -i°°  and  only  the  fourth  root  would  tend  to  -*-/».  Note  that  the  fourth  root  (let  us  designate  it  to 
be  kj  corresponds  to  the  waveguide  mode  which  propagates  in  the  opposite  direction  of  the  beam  velo¬ 
city.  If  the  movement  of  the  roots  is  in  accordance  with  that  depicted  in  Fig.  (6b).  a  convective  insta¬ 
bility  sets  in  when  c  is  slightly  increased.  On  the  other  hand,  if  the  roots  move  according  to  Fig.  (6c). 
an  absolute  instability  would  set  in  when  «  is  raised  to  a  higher  level.  In  other  words,  an  absolute  insta¬ 
bility  sets  in  only  if  the  merging  of  roots  in  the  k  plane  involves  the  backward  wave  mode  of  the 
waveguide.  This  result  perhaps  is  anticipated  from  the  outset  because  the  internal  feedback,  which  is 
required  for  an  oscillation,  can  be  provided  for  only  by  the  wave  supported  by  the  wavegutu.  which 
propagates  in  the  reverse  direction  of  the  beam.  This  wave  can  be  excited  only  if  the  beam  strength  is 
so  high  that  all  modes  are  strongly  coupled. 


•We  assume  that  the  solutions  assume  a  dependence  of  esrW  -  /*:). 
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From  these  considerations,  we  then  see  that  the  merging  of  k  roots  at  k  —  8 bc  in  Fig.  '5b»  must 
correspond  to  the  outset  of  a  convective  instability  since  the  merging  roots  arise  from  the  beam  mode 
and  the  forward  wave  mode  of  the  waveguide.  This  convective  instability  sets  in  as  long  as  c  >  0.  in 
agreement  with  the  conclusion  of  Appendix  A.  As  c  is  raised  slightly  above  zero,  we  encounter 
another  double  root  of  k  [cf.  kA  in  Fig.  (5c. d)].  This  double  root  of  k.  however,  does  not  correspond  to 
an  absolute  instability  since  the  merging  of  the  roots  does  not  involve  the  backward  wave  of  the 
waveguide  mode  either. 


Examination  of  Fig.  5  shows  that  the  only  multiple  root  which  involves  both  the  backward  wave 
and  the  (forward)  beam  mode  is  as  shown  in  Fig.  (5f).  Thus  the  onset  of  absolute  instability 
occurs  when  t  -  tCf,  where  tcf  is  given  by  [cf.  (Cl)] 


-  P?D(wsrkst) 

and  <i)  „ . k,t  are  the  real  root  determined  from  the  conditions 

-  o. 


bD  1 
bk 


b2D 

b:k 


o. 


A  little  algebra  yields  the  relevant  solution 


(C2) 


(C3) 

(C4) 


b0  +  6>/2fl,r 

<"1*  "  1  +  8,8;?  ’ 

k,t  -  ~  (wss  -  bp). 

These  equations  are  just  in  Eqs.  (9).  (8)  of  the  main  text.  Substitution 
gives  Eq.  (7). 


(C5) 

(C6) 

of  <C5)  and  (C6)  into  (C2) 


A  few  remarks  concerning  our  calculations  are  in  order.  It  can  readily  be  shown  froth  Eq.  (9)  that 
is  always  less  than  unity.  This  fact  that  oscillation  may  occur  at  a  frequency  lower  than  the  cutoff 
frequency  of  the  waveguide  is  hardly  a  surprising  result.  It  is  weil-known":>  that  when  a  waveguide  is 
loaded  with  a  plasma,  modes  with  a  frequency  below  the  cutoff  frequency  of  the  waveguide  can  exist. 
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In  the  present  case,  the  waveguide  is  partially  loaded  with  the  beam,  which  itself  constitutes  an  active 
medium.  This  brings  out  another  point  which  at  firs;  sight  may  also  seem  a  bit  paradoxical.  At  the 
grazing  condition,  one  may  then  ask.  how  can  an  absolute  instability  ever  exist  as  the  beam  mode  and 
the  waveguide  mode  both  propagate  in  the  forward  direction?  The  answer  is  that  at  very  low  beam 
current,  both  modes  are  weakly  coupled  and  indeed  no  absolute  instability  exists  according  to  the  usual 
notion  of  coupling  of  modes.  However,  as  a  sufficiently  high  current  level,  both  modes  are  strongly 
coupled.  In  fact,  all  modes  are  coupled,  including  the  waveguide  mode  which  propagates  in  the  oppo¬ 
site  direction  of  the  beam.  The  excitation  of  the  latter  modes  provides  an  internal  feedback  and  is  the 
basic  reason  for  the  possible  excitation  of  an  absolute  instability.  [Cf.  Fig.  ( 5f )  ] . 
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Tabie  I.  The  dimensioniess  frequency  (12  !r  and 
wavenumber  \k„  i  a-,  the  onset  of  absolute 
instability  t«  -  t.ci  for  various  3  Magnetic 
field  is  maintained  at  the  gracing  condition. 


8 

*>c 

<»!( 

Mr 

0 

1 

1 

0 

0 

.1 

0.995 

0.99985 

0.01216 

0.008'* 

0.9798 

0.9994 

0.0245 

0.0250 

0.9539 

0.9986 

0.03*2 

0.0465 

0.9165 

0.9974 

0.0506 

0.0*29 

■N 

0.8660 

0.9Q5S 

0.0649 

0.1046 

.6 

0.8000 

0.9935 

0.0806 

0.1423  ■ 

■* 

0T141 

0.9902 

0.0986 

0.1881 

.8 

0.6000 

0.9854 

0.1204 

0.2455 

0.4359 

0.9771 

0.1503 

0.3251  i 

Table  II.  Data  for  the  transition  to  absolute  instability 
at  various  magnetic  field  strength  (b,j  bp)  with 
&  fixed  at  0.266.  Shown  are  the  normalized  characteristic 
frequency  (ws).  wavenumber  (k,), 
threshold  current  (I)  and  the  axial  wavelength 
(X.)  at  the  transition. 


b'"b  o 

WS 

I  (amps) 

X.  (cm)  i 

.92 

.990471 

.097383 

45.42 

8.98  ; 

.94 

.99327 

.081894 

22.71 

10.68  i 

.96 

.995634 

.065996 

9.58 

13.25  i 

!  .98 

.997531 

.049659 

3.07 

17.613  ' 

|  1.00 

.998920 

.032845 

0.588 

26.625  - 

1  1.02 

.999759 

.015514 

0.02925 

56.377  | 

Table  III.  Data  for  transition  to  absolute  instability 
at  various  values  of  wall  resistivity  (8/r, ). 

Here  0,  -  .266.  b/b0  —  .98. 

Note  that  if  6  ^  0.  ks  is  complex,  and  that  <2,  <  1. 


6/r. 

ks 

I  (amps) 

0 

.997531 

(.049659.  0) 

3.07 

.002617 

.996291 

(.036153.  -  .05302) 

5.567 

.005234 

.996123 

(.02889,  -  .067291) 

7.684 

.007851 

.996193 

(.023109.  -  .07763) 

9.955 

.010468 

.996391 

(.018168,  -  .08602) 

12.391 

A- 2  3 


NRL  N^EMORANDL'M  REPORT  4303 


LU 

a 

O 


A-  2  5 


(  ig  2  —  The  waveguide  imute  ami  the  beam  mode. 
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Threshold  curreni  (|)  as  a  function  of  skin  depth  (wall  resistivity)  at  various  magnetic  field  strength  <  h  /  ) 

Here  fin  -  0.266,  <u,./2w  -  34.3CH:  and  \  -  0.533  cm. 
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(b) 
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(c) 


Trajectories  of  the  roots  in  the  k  plane  as  <u  changes  from  »,  to  »,-/«.  In  (b).  a  convective  instability  exists 
In  (c)  k,  corresponds  to  the  onset  of  an  absolute  instability. 
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EFFECTS  OF  VELOCITY  SPREAD  AND  WALL  RESISTIVITY 
ON  THE  GAIN  AND  BANDWIDTH  OF  THE 
GYROTRON  TRAVELLING-WAVE  AMPLIFIER 

I.  INTRODUCTION 

The  onset  of  oscillations  at  a  relatively  low  current  level  remains  a  major  obstacle  in  the  develop¬ 
ment  of  high  power  gyrotron  travelling  wave  amplifier  (gyro-TWA)11-2'  Analogous  to  the  conventional 
travelling  wave  tube,  it  has  been  experimentally  shown  that  the  introduction  of  wall  loss,  either  in  the 
form  of  a  lumped  or  distributed  loss'31,  may  render  the  gyro-TWA  stable  against  oscillations  resulted 
from  reflections.  Parallel  theoretical  studies14'  also  reveal  that  a  small  wall  resistivity  may  reduce  the 
danger  of  exciting  an  absolute  instability, (5>  whose  presence  can  lead  to  oscillations  even  in  the  absence 
of  reflections.  While  wall  resistivity  may  stabilize  the  amplifier  against  both  types  of  oscillations,  it  may 
also  modify  the  gain  and  the  bandwidth.  On  the  other  hand,  there  is  always  the  interesting  possibility 
that  it  may  also  lead  to  some  type  of  resistive  instability  which  would  not  exist  if  there  is  no  wall  resis¬ 
tivity.  Thus,  in  this  work,  we  study  in  detail  the  effects  of  a  distributed  wall  resistivity  on  the  gain  and 
bandwidth  of  the  gyro-TWA,  with  the  velocity  spread  in  the  electron  beam  accounted  for.  Previ¬ 
ously16,  the  effect  of  velocity  spread  has  been  studied  in  the  absence  of  wall  resistivity.  Here,  we  study 
the  combined  effects  of  the  wall  resistivity  and  velocity  spread  with  emphasis  on  the  former. 

Under  the  realistic  assumptions  that  the  skin  depth  8  of  the  wall  material  is  much  less  than  the 
wall  thickness  and  the  wall  radius  rw,  we  generalize  the  dispersion  relation  derived  earlier17'  to  include 
the  effect  of  wall  resistivity  and  of  velocity  spread.  We  analyze  the  gain  curves  over  a  wide  range  of 
parameters.  The  principal  results  are  summarized  as  follows. 

The  presence  of  a  small  resistivity  (e.g„  5/ r%  -  0.0026,  as  in  the  NRL  experiments),  only  slightly 
reduces  the  gain  of  the  amplifier.  More  specifically,  the  reduction  in  gain  is  approximately  1/3  of  the 
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cold-tube  loss  for  frequencies  in  the  neighborhood  of  the  grazing  intersection.  This  resuit  may  be  para¬ 
phrased  as  that  the  resistive  wall  is  3  times  as  effective  in  damping  the  backward  waves  reflected  from 
the  end  of  the  waveguide  than  the  forward  amplifying  waves,  leading  to  the  distinct  possibility  that  wall 
resistivity  may  stabilize  the  amplifier  without  significantly  reducing  its  gain.  (Recall  in  this  regard  that 
the  resistive  wall  may  also  stabilize  the  absolute  instability,  as  mentioned  earlier). 

Of  equal  interest  is  our  finding  that  wall  resistivity  can  increase  the  bandwidth  for  realistic  values 
of  &/>„..  This  increase  in  bandwidth,  which  is  quite  substantial  when  the  applied  magnetic  field  is  main¬ 
tained  exactly  at  the  grazing  value,  results  primarily  from  the  leveling  of  the  gain  curve.  The 
bandwidth  is  less  sensitive  to  wall  resistivity  if  the  magnetic  field  is  below  the  grazing  value.  The  pres¬ 
ence  of  wall  resistivity  also  broadens  the  unstable  band  in  the  frequency  spectrum:  growth  of  waves 
below  the  waveguide  cutoff  frequency  is  found  possible,  leading  to  some  optimism  that  very  wide  band 
amplifier  may  be  realized. 

A  modest  amount  of  velocity  spread,  say,  on  the  order  of  10  per  cent,  does  not  reduce  the  peak 
gain  significantly.  However,  the  bandwidth  may  be  significantly  reduced.  The  reduction  in  bandwidth 
is  most  noticeable  when  the  magnetic  field  is  exactly  at  the  grazing  value.  For  a  slightly  lower  value  of 
magnetic  field,  the  reduction  in  bandwidth  is  minimal,  even  at  a  current  level  as  high  as  10  amps. 

The  extensive  numerical  study  of  gain  and  bandwidth  points  to  the  following  optimum  design 
parameters:  A  beam  current  on  the  order  of  10  amps,  an  applied  magnetic  field  about  2%  below  the 
grazing  value,  a  wall  resistivity  with  a  skin  depth  about  five  thousandth  of  the  wall  radius.  For  a  gyro- 
TWA  with  such  design  parameters,  a  gain  on  the  order  of  3.5  db/wavelength  with  a  bandwidth  of 
approximately  4  per  cent  could  be  achieved.  Such  an  amplifier  may  operate  stably  with  significant 
efficiency. 

In  keeping  with  the  operation  of  the  amplifier,  namely,  amplification  of  a  signal  of  frequency  <u.  in 
this  paper,  we  shall  calculate  the  complex  axial  wavenumber  k.  as  a  function  of  (real)  w  from  the 
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dispersion  relation.  We  should  remark  that  in  so  doing,  care  must  be  exercised  to  obtain  and  to  inter¬ 
pret  the  gain  curve,  especially  when  the  current  exceeds  tne  threshold  value  for  the  onset  of  absolute 
instabilities. 141  The  proper  branch  of  the  solution  must  be  chosen  [cf.  Ref.  (5).  p.  38].  On  the  other 
hand,  calculations  of  complex  for  real  values  of  k.  might  easily  be  misused  in  connection  with 
amplifiers  (cf.  Ref.  (5),  p.  167], 

In  Section  II,  we  shall  formuiaie  the  dispersion  relation  including  the  effect  of  wall  resistivity  and 

of  the  velocity  spread  in  the  electron  beam.  The  numerical  results  are  presented  in  Section  III.  The 
major  results  and  implication  are  discussed  in  Section  IV. 

II.  FORMULATION 

The  inclusion  of  a  small,  but  finite  wall  resistivity  in  the  waveguide  requires  some  modification  of 
the  derivation  of  the  dispersion  relationship.  Considerable  simplification  results  when  two  assumptions 
are  introduced:  (1)  The  beam  is  sufficiently  tenuous  that  the  wave  with  which  it  interacts  is  primarily 
the  vacuum  field  supported  by  the  (lossy)  waveguide;  (2)  The  thickness,  as  well  as  the  inner  radius  r„ 
of  the  resistive  wall  are  much  larger  than  the  skin  depth  8  associated  with  the  wall  material. 

The  first  assumption  was  adopted  in  virtually  all  literature  on  gyrotron.  The  excellent  agreement 
between  experiments  and  the  theory  developed  under  this  assumption  leaves  little  doubt  regarding  its 
validity.  Implied  also  in  this  assumption  is  that  the  beam  mode  and  the  waveguide  mode  are  close  to 
synchronism.  One  should  take  precaution,  therefore,  if  the  present  formalism  is  to  be  extrapolated  to 
explore  the  possibility  of  designing  an  extremely  wideband  amplifier,  even  though  such  a  possibility  is 
potentially  attractive.  The  second  assumption  is  also  generally  valid  for  any  realistic  gyro-TWA  experi¬ 
ment.  For  example,  in  the  planned  NRL  experiments  r»  -  .533  cm,  8  =*  1.4  x  10~3  cm  at  35  GHz.(3) 
Thus  8 / rw  a  2.6  x  10-3,  an  extremely  small  quantity  compared  with  unity. 

Our  calculation  of  the  dispersion  relationship  follows  closely  with  Chu.  Drobot,  Szu.  and  Spran- 
gle.,7>  First,  the  vacuum  fields  of  the  TE<,„  mode  in  a  cylindrical  lossy  waveguide  are 


if  8  -  0 
if  6  *  0 


To  make  full  use  of  the  mathematical  property  of  Bessel  functions,  let  us  define 

_5_  (1  -  j) 

r„  2 

Thus 


J l  < *>  “  0. 


(8) 

(9) 

(10) 

(11) 
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For  ready  comparison  with  the  cold  tube  loss  rate  resulting  from  wall  resistivity,  we  calculate  the 
axiai  wavenumber  k.  for  the  emply  waveguide,  From  Eqs.  (5)  and  to),  we  find,  for  the  empty 
waveguide,  k.  -  k._£  where 

-  <4  -  a-„  (1  -  /)  — .  (12) 

wuh 

li/: 

i.  —  OJ*  ; 

r^20  *  &  on  ( ) 

f* 

being  the  axial  wavenumber  of  the  corresponding  lossless  waveguide.  Note  that  Eq.  (12)  is  valid  even 
if  Lc  -  0.  It  will  be  used  to  calculate  the  attenuation  rate  of  the  waves  in  the  cold  tube  due  to  the  wall 
loss.  On  the  other  hand,  measurements  of  attenuation  rate  for  the  cold  tube  at  various  frequencies  may 
be  used  to  determine  the  skin  depth  of  the  wall  material  by  Eq.  (12). 


We  now  consider  the  interaction  of  the  electron  beam  with  the  above  electromagnetic  waves.  In 

the  typical  configuration  of  a  gyro-TWA  (Fig.  1),  the  electrons,  forming  an  annular  beam,  are  guided 
by  a  uniform  magnetic  field  ( B0  e.).  They  move  along  helical  trajectories.  Ideally,  all  electrons  have 
the  same  perpendicular  velocity  vi0  and  the  same  parallel  velocity  v.„ ,  with  their  guiding  centers  uni¬ 
formly  distributed  on  a  surface  of  constant  radius  r0  (Fig.  1 ) .  The  electromagnetic  fields  represented  by 
Eqs.  ( 1 ) - (3)  induce  a  current  perturbation  in  the  beam.  This  response  in  the  current  has  been  com¬ 
puted  by  Chu,8)  for  the  lossless  waveguide  by  solving  the  linearized  Vlasov  equation.  His  calculation 
can  readily  be  extended  to  the  present  case  of  a  lossy  waveguide  because  the  field  solutions  y0(a„r)  and 
J\(anr)  in  Eqs.  (l)-(3)  are  entire  functions  of  a„r.  In  order  that  this  perturbation  current  Jtv  be  con- 
sistant  with  the  field  solutions  (l)-(3),  they  must  be  related  by  the  Maxwell  equations  as  follows: 


where  4"’  is  the  0  component  of  the  perturbation  current,  Strictly  speaking,  Eq.  (14)  is  an  integral- 
differential  equation  since  4<n  is  a  functio  n  1  of  the  fields  which  are  to  be  solved  self-consistentlv. 
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The  first  assumption  introduced  at  the  beginning  of  this  section  can  now  be  used  to  simplify  the 
formulation.  By  mis  assumption  B.‘  is  given  by  Ec  <1;  Multiply  equation  « 14 1  by  r/-.(Q,r'  and 
integrate  the  resultant  equation  from  r  -  0  to  r.  *  We  men  obtain 


in' 

— : —  a:.*  —  a  p 
c* 

where  we  have  used  equation  ( 1 1  >  andar.  \  *  are  given  by  equations  (7)  and  (10)  respectively.  This 
procedure  is  justified  as  all  quantities  in  equation  <  1 4 >  are  entire  functions  of  r.  For  the  same  reason, 
the  modified  dispersion  relation  now  reads 


8rra .  e 


~  ik  . .  'iui; 


c  r.  **  J ;  i;„ 


dr  r  Jl'1'  J j  ( a.r ) 


(15  i 


where  ft  c  - 


OT 

c: 


-  a  *  -  TT^rr-r  J0  p.  {p-  ■*  > 

4  on*' (j  '4  on ' 

(to :  —  k}c:)p‘Hs(arril.anrL  ) 

•yJ/n*c:(<o  —  fcv.  —  sft.)* 

(to  -  k.\.)Qs[a rr0.a„rL  ) 

■y(to  —  A.v.  —  sfl,.) 


(16) 


eB0/ymc,  rL  -  v/ft..  v  s  .V>:/mc:  is  the  dimensionless  beam  density  parameter.  .Vis  the 


number  of  electrons  per  axial  length,  and  the  functions  //,  and  Qs  are  defined  by 


H%  (x,y )  =  [7,  (xl/'j  (.v )]: 
and 


&(*,>)  S  2HS (xv )  -r  yj't  0  )/", (y ){J‘ Or  H 1  -  s'/x') 

—  [y,U)]J)  -  2s:J,  (x)J\  (x  )J\  (y )  [y/',  (y ) 

-  Js(y)]/xy. 


In  Eq.  (16),  g(p,.p;)  is  an  arbitrary  function  of  the  momentum  p_  and  p.  satisfying  J  g(p..p;) 
d:'p  «  1.  The  most  suitable  electron  beams  for  the  gyrotrons  are  those  generated  from  a  magnetron- 
type  electron  gun.  Typically,  such  beams  are  characterized  by  a  negligible  energy  spread,  hence  the 
velocity  spread  comes  primarily  from  the  pitch  angle  spread  of  the  electron  velocity.  To  model  this 
kind  of  beam,  we  let 


g(Pi.P:) 


K  5(y— y„)  exp 


-(£•  ~ 
(AP;)J 
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where  A'  is  a  normalization  constant  and  £p:  is  the  mean  axtai  velocity  spread,  and  6  is  the  delta  func¬ 
tion  of  Dirac. 

Equation  (16*  is  identical  in  form  with  that  for  the  lossless  waveguide"-  6'  except  that  a  r  now  has 
a  small  imaginary  pan.  and  that  we  have  neglected  the  azimuthal  modes.  It  incorporated  two  obvious 
physical  effects  of  the  wall  resistivity:  (a!  its  modification  of  the  waveguide  mode  and  hence  (b)  the 
modified  beam-wave  interaction  as  a  result  of  the  modification  in  the  electromagnetic  field.  The  firs: 
effect  is  represented  by  the  complex  a„  appearing  in  the  LHS  of  Eq.  (16)  and  the  second  effect  by  that 
in  the  RHS.  In  Appendix  B.  we  show  that  the  second  effect  is  negligible  in  comparison  with  the  first 
effect.  Thus  we  may  replace  a„  by  aBn  in  the  right  hand  side  of  Eq.  (16). 

It  is  anticipated  that  wall  resistivity  should  reduce  the  forward  gain  of  the  amplifier.  Actually,  the 
reduction  in  the  gain  is  considerably  less  than  that  expected  of  the  vacuum  waveguide.  In  Appendix  B. 
we  estimate  that  the  reduction  in  the  gain  by  wall  resistivity  is  approximately  1/3  of  the  cold  tube  loss 
(by  resistivity)  at  the  same  frequency.  This  analytic  result  is  supported  by  our  numerical  results.  It 
may  be  regarded  significant  since  it  means  that  wall  resistivity  may  stabilize  the  amplifier,  by  damping 
out  any  reflected  signal  from  the  end  of  the  waveguide,  without  sacrificing  the  forward  gain  as 

significantly. 

In  Section  III,  we  shall  present  numerical  results  for  various  combinations  of  6 />»,  velocity- 
spread.  beam  current  and  the  applied  magnetic  field.  The  results  are  obtained  from  the  full  dispersion 
relationship  (16),  and  from  its  modified  form  when  velocity  spread  is  present. 

III.  NUMERICAL  EXAMPLES 

The  general  dispersion  relationship  (16)  is  applicable  to  all  cyclotron  harmonics  (s)  and  to  all 
radial  mode  numbers  (n).  Thus,  a  huge  parameter  space  needs  to  be  explored  for  an  exhaustive  inves¬ 
tigation  of  the  performance  of  the  gyro-TWA.  To  limit  our  studies  of  the  effects  of  wall  resistivity  and 
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of  velocity  spread,  we  shall  only  examine  the  lowest  cyclotron  harmonic  of  the  TEo;  mode.  Ever.  then, 
four  parameters  need  to  be  specified  to  obtain  the  gain  curve.  They  are: 
ia>  the  magnetic  field  B0,  in  units  of  the  grazing  value  Bq. 
tb>  the  current /. 

tc)  the  wall  resistivity,  measured  by  6,  r*.  and 
<d)  the  velocity  spread,  measured  by  Ip.!  p.0. 

The  magnetic  field  is  said  to  assume  the  grazing  value  when  the  non-relativisitic  electron  cyclotron  fre¬ 
quency  equals  to  the  cut-off  frequency  of  the  waveguide.  The  numerical  examples  obtained  may  then 
be  used,  after  properly  scaled,  to  predict  the  performance  of  the  amplifier  operated  at  higher  harmonics 
and  at  higher  order  modes. 

To  facilitate  comparison  of  results,  at  the  expense  of  some  duplications,  we  shall  present  our 
numerical  calculations  in  three  subsections.  In  each  subsection,  gain  curves  are  given  in  which  three  of 
the  four  parameters  mentioned  above  are  held  fixed  and  only  the  remaining  one  is  allowed  to  vary.  All 
gain  curves  are  drawn  with  the  same  scale.  The  individual  effects  are  thus  placed  in  sharp  focus.  In 
Subsection  a,  in  each  graph  of  the  gain  curves,  only  the  wall  resistivity  is  varied.  In  Subsection  b.  only 
the  velocity  spread  is  varied  and  in  Subsection  c,  only  the  current  is  varied.  In  all  cases  studied,  we 
assume  that  j3„  -  0.266,  vy  v,l0  -  1.5  and  that  the  beam  energy  is  70  keV.  These  assumptions  closely 
represent  the  NRL  experiments.11’  We  shall  postpone  to  the  next  section  for  some  general  discussion  of 
the  role  of  resistivity  in  instabilities. 

a.  The  Effect  of  Wall  Resistivity 

In  this  subsection,  we  present  the  gain  curves  in  which  the  wall  resistivity  (5/rJ  is  allowed  to 
vary.  We  shall  plot  the  normalized  spatial  growth  rate  k,  =  -ktrw  as  a  function  of  normalized  fre¬ 
quency  SI  s  w  rjc.  For  a  lossless  waveguide,  SI  assumes  the  value  3.832  at  the  cutoff  frequency  and 
3.975  at  the  grazing  frequency.  In  both  Figs.  (2)  and  (3),  we  set  B/Bq  -  1  and  lp;  “  0.  In  Fig.  (2). 
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I  —  1  amp  and  in  Fig.  3.  /  -  3  amps.  Also  shown  in  parentheses  in  these  figures  are  the  bandwidths. 
which  represent  the  frequency  range  inside  which  the  gain  is  within  85°c  of  the  peak  value.  It  is  seen 
from  these  two  figures  that  the  gam  is  not  reduced  by  a  significant  amount  as  5/r.  is  raised  from  0  to 
0.0025  (again,  in  the  NRL  experiments.  5/r„  =  0.0025).  If  6/ru  is  further  raised  from  0.0025  to  0.01. 

amounts  to  an  increase  of  wall  resistivity  by  a  factor  of  16.  we  see  that  the  bandwidth 
significantly  increases  [c.f.  Fig.  3]  while  the  gain  is  not  substantially  reduced.  In  fact,  if  we  compare  the 
numerical  value  of  this  reduction  in  gain  with  the  cold  tube  loss  as  determined  from  Eq.  (12)  [c.f.  Fig. 
41  the  former  is  approximately  1/3  of  the  latter,  over  the  useful  bandwidth  in  the  vicinity  of  the  grazing 
frequency.  This  result  is  in  qualitative  agreement  with  the  analytical  study  given  in  Appendix  B.  The 
sharp  peak  in  the  5 />,  *  0  curve  in  Fig.  3  is  indicative  of  the  presence  of  an  absolute  instability'4'  in 
that  case. 

In  Figs.  5-7,  we  include  the  effect  of  the  velocity  spread  while  the  magnetic  field  is  still  main¬ 
tained  at  the  grazing  level.  It  is  again  seen  that  the  bandwidth  increases  with  wall  resistivity,  and  the 
forward  gain  is  not  reduced  significantly. 

In  Fig.  8.  we  set  B/BG  —  0.98,  7-7  amps,  and  A p.  -  0.  The  bandwidth  also  increases  with  wall 
resistivity,  but  its  dependence  on  resistivity  is  less  sensitive  than  that  in  the  case  of  exact  grazing  mag¬ 
netic  field.  The  slight  decrease  in  bandwidth  as  8/r*  changes  from  0.005  to  0.01  as  shown  in  this  figure 
suggest  the  presence  of  resistive  damping. 

It  is  of  interest  to  point  out  that  even  though  the  peak  gain  is  reduced  by  wall  resistivity,  the 
unstable  band  of  the  amplifier  is  enlarged.  Virtually,  in  all  cases,  the  unstable  band  extends  to  frequen¬ 
cies  below  the  cut-off  frequency  of  the  waveguide  if  wall  resistivity  is  included,  thus  suggesting  the 
presence  of  some  resistive  wall-induced  instability.  We  shall  comment  on  this  dual  role  of  resistivity  in 
Section  IV. 
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b.  The  Effect  of  Velocity  Spread 

In  this  subsection,  we  evaluate  the  effect  of  velocity  spread  on  the  gain  and  bandwidth  of  the 
gyro- TWA.  In  each  figure,  all  parameters  are  held  fixed  except  the  velocity  spread  \p:i  pzc.  We  take 
the  values  of  lp:/pz0  to  be  0.  7%  and  15%.  It  is  found  that  the  velocity  spread  reduces  the  gain  and 
bandwidth  significantly  when  the  magnetic  field  is  maintained  at  the  grazing  level.  On  the  other  hand, 
when  Bi Be  <  1,  velocity  spread  does  not  seem  to  reduce  the  gain  nor  the  bandwidth  to  as  a  great 
extent,  even  at  a  high  current  level.  [Of  course,  one  might  also  worry  that  at  a  high  current  level,  the 
beam  quality  may  deteriorate  and  the  velocity  spread  may  increase]. 

In  Figs.  9-11.  we  set  B/Bc  —1.7—1  amp.  and  let  6/r.  —  0.  0.0025  and  0.1  respectively.  We 
see  that  at  this  low  current  level,  the  presence  of  velocity  spread  reduces  the  gain  and  bandwidth  quite 
significantly  [c.f.  Fig.  11].  This  remains  to  be  the  case  when  the  current  is  raised  to  3  amps,  as  shown 
in  Figs.  12-14,  where  we  still  set  B  »  Ba¬ 
li  we  lower  the  magnetic  field,  so  that  B/Bq  *  0.98,  [c.f.  Figs.  (15M17)],  then  the  velocity 
spread  has  a  lesser  influence  in  the  gain  and  in  the  bandwidth.  In  these  figures,  we  have  set  8/ /•»  - 
0.0025.  A  higher  value  of  8/r„,  such  as  0.01  as  used  in  Figs.  18-21.  would  permit  the  amplifier  to 
operate  stably  at  substantially  higher  current.141  It  is  apparent  from  these  figures  that  the  bandwidth,  as 
well  as  the  gain,  again  does  not  greatly  suffer  from  the  presence  of  velocity  spread. 

c.  The  Effect  of  Increasing  Current 

Finally,  we  examine  the  dependence  of  gain  and  bandwidth  on  the  current.  When  the  current  is 
sufficiently  high,  an  absolute  instability  sets  in.  There  would  then  be  some  ambiguity15'  regarding  the 
proper  solutions  of  the  amplifying  waves.  Such  an  ambignity  arises  from  the  intersection  of  the  branch 
cut  in  the  <u  plane  with  the  real  «  axis.  It  does  not  seem  to  pose  a  serious  problem  in  the  present  study 
when  the  beam  current  exceeds  the  threshold  current  only  by  a  moderate  factor,  such  as  1.5.  We  have 
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exercised  some  care  in  choosing  the  correct  solution,  and  we  may  refer  to  Ref.  U)  for  a  detailed  discus¬ 
sion  of  the  absolute  instability  3nd  for  the  calculation  of  the  threshold  current  in  the  gyro-TWA. 

Here,  we  report  that  whiie  the  gain  increases  with  the  current,  as  expected,  the  bandwidth 
decreases  once  the  current  exceeds  the  threshold  value  for  the  onset  of  absolute  instability.  This  may¬ 
be  intepreted  as  the  tendency  of  the  system  to  respond  mostly  to  a  particular  frequency,  namely,  the 
natural  frequency  of  oscillation  when  the  absolute  instability  is  present,  thereby  leading  to  a  narrow 
bandwidth  that  results  from  the  sharp  peak  in  the  gain  curve  in  the  neighborhood  of  that  natural  fre¬ 
quency.  These  statements  are  corroborated  in  Figs.  22  and  23.  where  the  gain  curves  tend  to  peak  at 
the  cutoff  frequency  (i.e.,  wc  -  <ucrjc  *  3.832),  which  also  is  the  frequency141  for  the  onset  of  the 
absolute  instability.  From  Ref.  (4)  the  threshold  current  for  parameters  shown  in  Fig.  22  is  1.5  amps 
and  that  for  Fig.  23  is  5.0  amps. 

On  the  other  hand,  if  the  current  is  below  the  threshold  value,  then  increasing  the  current  would 
lead  to  the  increase  of  both  gain  and  bandwidth,  as  clearly  indicated  in  Figs.  24  27.  In  these  figures,  we 
set  B!Bq  —  0.98  so  that  the  threshold  current  is  high*4’.  The  threshold  current  for  both  Figs.  24  and 
25  is  5.5  amps  and  that  for  both  Figs.  26  and  27  is  12  amps.  The  velocity  spread  in  Fig.  24  is  7%  and 
that  in  Fig.  25  is  15%.  Figure  26  and  Fig.  27  differ  also  only  in  velocity  spread.  The  curves  in  Fig.  26 
are  hardly  modified  if  the  velocity  spread  in  7%  instead  of  zero  [c.f.  Figs.  15-21], 

IV.  DISCUSSIONS 

One  thing  that  stands  out  among  the  numerical  results  is  the  dual  role  played  by  electrical  resis¬ 
tivity:  (a)  in  its  reduction  of  the  gain  of  the  cyclotron  maser  instability  and  (b)  in  its  excitation  of 
amplifying  waves  in  a  regime  where  such  amplification  would  have  been  impossible  if  the  resistivity  is 
absent.  Such  a  dual  role  of  dissipative  effects  is  of  common  occurrence  in  physical  systems  with  avail¬ 
able  free  energy.  If  in  such  a  system  an  instability  of  the  reactive  type15'  is  already  present  even  in  the 
absence  of  dissipative  effects,  the  introduction  of  dissipation  would  reduce  the  growth  rate  becausJ  it 
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takes  away  the  free  energy  that  drives  the  instability.  On  the  other  hand,  if  there  is  no  instability  m 
such  a  system  in  the  absence  of  dissipation,  the  introduction  of  dissipative  effects  would  lead  to  growth 
of  waves,  as  dissipation  would  induce  relaxation  of  the  system  which  is  on  the  verge  of  becoming 
unstable  to  begin  with.  Numerous  examples  of  these  phenomena  can  be  found  in  the  literature  of  clas¬ 
sical  mechanics,  hydrodynamics,  plasma  physics,  electron  beam  devices,  astrophysics,  etc. 

The  major  results  and  their  implication  of  the  present  study  may  then  be  summarized  as  follows: 

(1)  In  terms  of  power,  gain,  bandwidth,  efficiency,  and  stability,  the  gyro-TWA  would  achieve 
best  overall  performance  when  the  magnetic  field  is  maintained  just  slightly  (say.  2%) 
below  the  grazing  value.  Stable  operation  with  significant  gain  and  bandwidth  (such  as  3-4 
db  per  wavelength  in  gain  and  4%  in  bandwidth)  may  be  achieved  at  a  current  of  10  amps, 
or  even  higher.  if5/r»  ^  0.005. 

(2)  Wall  resistivity  leads  to  a  reduction  in  gain  by  an  amount  on  the  order  of  1/3  of  the  cold 
tube  loss.  It  leads  to  an  increase  in  bandwidth  when  the  magnetic  field  is  maintained  at 
the  grazing  value.  If  the  magnetic  field  is  below'  the  grazing  value,  and  if  the  current  is 
low,  the  presence  of  wall  resistivity  may  actually  reduce  the  bandwidth.  The  latter  may  be 
regarded  as  an  example  of  resistive  damping. 

(3)  Velocity  spread  typical  of  present  electron  guns  (Ap-/p-0  =  5-  10%)  does  not  significantly 
reduce  the  peak  gain.  The  underlying  reason  is  that  its  effect  enters  in  the  combination 
k.-(Ap.)2  qualitatively  and  the  amplifier  reaches  a  peak  gain  when  w  =  tuc,  at  which  k.  is 
small.  At  higher  frequencies,  k.  becomes  sizable  and  the  effect  of  velocity  spread  would 
then  be  important.  Thus,  the  bandwidth  may  be  significantly  reduced,  especially  if  the 
magnetic  field  is  maintained  at  the  grazing  value. 

(4)  The  gain  and  bandwidth  of  the  amplifier  are  less  sensitive  to  velocity  spread  or  to  resis¬ 
tivity  when  the  magnetic  field  is  below  the  grazing  value.  Thus,  for  stability  of  the 
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amplifier,  we  need  a  wail  resistivity  whose  skin  depth  is  on  the  order  of  5  x  1C r-  r. .  To 
reduce  the  deteriorating  influence  of  the  velocity  spread  on  the  bandwidth,  the  magnetic 
field  should  be  maintained  slightly  below  the  grazing  value.  The  current  should  be  on  the 
order  of  10  amps,  which  is  sufficiently  high  for  gain  and  power,  yet  low  enough  to  avoid 
oscillations.  These  arguments  form  the  basis  of  the  parameters  quoted  in  the  Abstract. 
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is  much  less  than  the  free  space  wavelength  c/u.  Using  (A4)  into  (Al).  and  introducing  the  transfor- 


S.'"-  r~l  -  <t> . 


we  rewrite  equation  (Al)  for  Region  II  as 


d:d>  .  f.  .  1  1  ,  n 

— - +  — -  e>  -  0. 

dr *  |_  4r- 


<  r  <  oo . 


The  last  term  1/ 4r-  is  clearly  small  in  comparison  with  uiy.ccr  for  all  r  >  r,  and  hence  is  neglected.  In 
order  that  the  solution  to  (A7)  is  bounded  as  r  —  »,  we  have 

B.'u(r)  -  r-'! '■  6(r)  =  B  r"1'2  e"""’"6;  r  >  rR  (A8) 

where  B  is  an  arbitrary  constant.  Two  conditions  are  now  imposed  on  the  solutions  (A3)  and  (A8)  at 

the  boundary  r  —  r„  separating  the  two  regions: 


*"'(') I  -  -  B‘"(r)L 


dBl'Hr)  d3.m(r ) 

- -  -  -  (A10) 

Condition  (A9)  assures  the  continuity  of  the  magnetic  field  and  (A10)  the  continuity  of  the  tangential 
electric  field.  Observe  from  Eq.  (A8)  that  the  variation  of  the  magnetic  field  in  the  wall  region 
(r  >  rw)  is  dominated  by  the  exponential  factor  since  8  «  r„.  Thus,  in  imposing  condition  (A10), 
we  can  ignore  the  geometrical  factor  r~i;1  in  Eq.  (A8),  It  can  then  immediately  be  shown  that  (A9> 
and  (A10)  lead  to  the  dispersion  relationship 


6  1  -  / 


(*1,1 0  J0lknXrw)  -  /,(*„, rj 


which  is  reproduced  as  Eq.  (4)  in  the  main  text.  Of  course,  this  dispersion  relation  can  also  be 
obtained  by  first  formulating  the  field  solutions  in  the  wall  region  in  terms  of  Bessel  functions  (with 
complex  arguments),  to  be  followed  by  the  use  of  asymptotic  formulas  of  the  Bessel  functions  when 
S/r„  «  1.  The  formulation  presented  above  illustrates  clearly  the  nature  of  such  approximations. 
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With  ;  H  <r,,  r,  anc  5 >>.  «  1.  Eq.  (All)  aamus  the  approximate  solution 


where  J{  l£0„)  -  0  and  A;  is  given  by 


i  “  £  ar  +  -±i 


( A 1 2 ) 


(5/r.) 

1  -  i 

*7 

£  on 

(ay,(0/ac);-;M 


( A 1 3 ) 


Equation  (A12)  reduces  to  Eq.  (6)  of  the  main  text  after  we  use  in  (A13)  the  formula 
d-M{)/3£  -  J0{O  at  the  roots  lor  of  Jy 


We  summarize  by  stating  that  the  vacuum  fields  of  the  TE0„  mode  in  a  lossy  waveguide  are  given 
by  Eqs.  (l)-(3)  in  the  main  text.  The  radial  wavenumber  an  in  these  expressions  is  given  by  Eq.  (7). 
The  influence  of  the  wall  resistivity  on  the  beam  dynamics  will  be  discussed  in  Appendix  B. 
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EFFECTS  OF  RESISTIVE  WALL  ON  THE  DYNAMICS  OF 
THE  ELECTRON  BEAM 


In  this  Appendix,  we  assess  the  influence  on  the  electron  beam  dynamics  by  a  small,  but  finite, 
wall  resistivity.  We  shall  show  (a)  that  the  loss  in  the  forward  gain  of  the  amplifier  is  approximately 
1/3  of  the  cold  tube  loss  and  (b)  that  the  electron  beam  essentially  experiences  the  vacuum  field  inside 
a  lossless  waveguide.  Again,  our  estimates  will  be  based  on  the  smallness  of  8/  r „  and  we  shall  consider 
complex  roots  of  k.  for  real  values  of  w. 


Near  synchronism.  -  k?c‘  -  a~c-  as  0  and  <o  -  k. v.0  -  s£),  =  0.  The  second  term  in  the 
right  hand  member  of  Eq.  (16)  can  be  ignored  compared  with  the  first  term.  We  can  further  approxi¬ 
mate  w:  -  k?c'  by  a;c‘  in  the  right  hand  side  of  (16).  Thus,  Eq.  (16)  is  simplified  to  read 


Ur  -  ks  c:  —  a2c2)  (w  —  k.v:c  -  sflr): 
where,  for  analytical  tractability,  we  have  set  Ap.  -  0. 


4 

y„  ion  ^6  (£  on  ) 


Hs(a„r0\  a„rL) 


(Bl) 


With  5  =  to rj c,  k  —  k.rw ,  Eq.  (Bl)  becomes 


k2  -  Z2  +  ci„ 


1  +  —  (I  -  /) 


-  to  —  Q , 

k  - 


-£  =  £,+•  /£, 


(B2) 


where  we  have  used  Eq.  (7)  and  fif  S  iflrr»/c.  In  Eq.  (B2).  the  parameter  £.  whose  real  and  ima¬ 
ginary  parts  are  denoted  E,  and  £,,  is  proportional  to  the  right  hand  member  of  Eq.  (Bl).  It  is  clear 
from  (7)  and  (Bl)  that 


E, 


(B3) 


Thus,  both  sides  of  (B2)  are  modified  by  wall  resistivity  through  the  introduction  of  a  term  of  order 
5 /r„.  The  statement  (b)  at  the  beginning  of  this  Appendix  claims  that  £,  can  be  neglected.  The  argu¬ 
ment  supporting  this  claim  will  be  given  later  in  the  Appendix.  We  shall  now  prove  the  statement  (a). 
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In  the  absence  of  the  beam.  £ —  0.  and  the  wave  inside  the  lossy  waveguides  attenuates  at  a  rate 
given  by  the  imaginary  part  of  kw.  where 


k,  “  Ul'  -  i‘n  I  —  —  (1  —  ;)  =r  HZ'  —  1  ~ 


Cor  (/  -  1 )  6/r. 

-  ci„) 


The  last  expression  is  valid  once  to/ Cor  ~  1  >  6/r»„  i.e..  once  the  frequency  a,  exceeds  the  cut-off  fre- 
quencvtu.  by  the  (small)  amount  (6 //■„.) cof.  Thus  the  normalized  attenuation  rate  of  the  cold  tube  is 


Ifn  ( kw )  “  i  or  ^  ^ 


(=:-ciB)1;:. 


To  estimate  the  reduction  of  the  gain  by  the  wall  resistivity,  we  neglect  £  in  (B2)  and  rewrite 


(B2)  as 


where 


(k  -  kj  (k  -  kB)2' 


(k  4-  kj 


r  s  -  nf 

which  is  real  for  real  values  of  w.  If  we  consider  the  vicinity  of  the  grazing  condition,  at  which 

*.  ■ 

P II 

then  we  may  rewrite  Eq.  (B6)  as 


where  we  have  defined 


k'  =  k  -  kB 


and  kg  =  Z0„p„  (1  —  1/2  is  the  characteristic  wavenumber  at  the  grazing  condition.  Note  that 

lm(£')  -  Im(£).  Upon  using  (B8)  and  (B4) ,  we  have 


{on  I  8  I  (/*"“!) 


Equation  (BID  may  now  be  used  in  Eq.  (B9)  to  estimate  the  effect  of  a  small  (8/r.)  on  the  gain  per 
unit  length,  which  is  measured  by  the  imaginary  part  of  k'.  We  represent  the  solutions  to  Eq.  (B9)  as 
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where  k\.  is  the  solution  to  (Bqi  when  wall  resistivity  is  absent,  ti  e..  6  =  0).  Thus,  k  £.  satisfies 

k'c  -  E j  (B13) 

and  the  imaginary  part  of  k':  represents  the  reduction  of  the  forward  gain  by  the  wall  loss.  By  substi¬ 
tuting  (B12)  into  (B9)  and  by  assuming  k\  «  .'c'o,  we  obtain 


7.  -  1  ,7  7  ,  ion  <5/0  0-1) 

*  i  “  ~y  '■*b~  - _  - - - - 


(B14I 


Vw*  -  So'n 

whose  imaginary  part  is  approximately  1/3  of  lm(£j  by  Eq.  (B4).  These  qualitative  estimates  are 
confirmed  by  detailed  numerical  calculations  (c.f.  Fig.  (4)).  Thus,  the  reduction  in  the  forward  gain 
because  of  the  wall  resistivity  is  only  one  third  of  the  cold  tube  loss. 


We  finally  show  that  a  small  imaginary  part  in  the  RHS  of  the  dispersion  relationship  (B2)  [or 
(B9)]  has  a  negligible  influence  on  the  gain.  To  this  end.  let  us  consider  Eq.  (B9)  and  add  a  small  ima¬ 
ginary'  part  to  its  right  hand  side.  Then  (B9)  becomes 

[£'-‘•81)-  k'H\  +  /6:)  (BIS) 

where  5]  and  8:  are  both  of  the  order  of  6/r,  and  k'l  “  £].  Again,  expressing  the  approximate  solu¬ 
tion  to  (BIS)  in  the  form  of  (B12).  we  find 

=  — J  [5i  -  *’08;]-  (B16) 

Since  6;  and  6;  are  of  the  same  order  and  k'0  «*  £7'3  is  a  small  quantity  whose  order  of  magnitude  is 
typically  in  the  range  between  0.01  and  .1.  we  conclude  that  the  introduction  of  a  small  imaginary  part 
in  the  RHS  of  the  dispersion  relation  (16)  has  a  much  smaller  effect  than  that  appearing  on  the  LSH  of 
(16).  In  physical  terms,  the  modified  fields  by  the  wall  resistivity  does  not  significantly  affect  the  beam 
current  perturbation  in  response  to  these  fields. 
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M  B/5g  =  0.98 

I  =  14  amps 
,  5/rw  =  0.01 


Fig.  22  —  Normalized  amplification  rate  kt  versus  normalized  frequency  S.  Here, 
only  the  current  /  is  allowed  to  vary,  all  other  parameters  B/ Bc,  i/rw.  are 

fixed  at  the  values  as  indicated  in  the  figure  The  numbers  in  paranthesis  denote  the 
bandwidths 
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MODE  COMPETITION  IN 
GYRCTRON  CAVITY 
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ASSUMPTIONS 

Initially,  all  the  electrons  have  the  same 
and  v  (i.e.  cold  beam) 

Initially,  all  electrons  have  the  same  guiding 
center  position  rQ  (i.e.  thin  annular  beam*' 

Mode  frequency  and  spatial  structure  are 
unperturbed  by  the  presence  of  the  beam  (i.e. 
tenuous  beam) 

There  is  a  steady  state  in  the  presence  of  many 
modes  in  the  cavity. 

The  equations  of  motion  of  each  electron  are 
independent  of  the  other  electrons  (i.e.  no 
space  charge  effects) 
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EQUATIONS  OF  MOTION  FOR  EACH  ELECTRON 


dP 

dt 


El'  '!o  t  £  Si1 

i=l 


d 


dt 


e 

- «-  v 

mc‘  — 


/  i  =  1,  2, 


M 


where 


\  =  £l  +  P2/m:c2J  2 
dx 

v  =  =  — 


Initial  values  of  A'v  .  (t  )  =  0. 

i  o 

For  the  mode  TEon^,  the  EM  field  is  equal  to: 

EQ  =  E  J,  (k  r)  sink  z  cos(u)t  +  4  ) 

0  o  l  n  z  o 

kz 

B  =  —  E  J,  (k  r)  cosk  z  sin(oot  +  <t) 
r  u;  o  1  n  z  o 

k 

B  =  -—  EJ  (kr)  sink  z  sin(u>t  +  $  ) 
z  oo  o  o  n  z  o 


where 


J,  (k  r  )  =0 
1  n  w' 


k  = 
z 


TT  £ 


,  £  =  1,2, 


w 


'  [kn  +  KT  c 
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MAIN  QUANTITIES  EVALUATED 

_K .  (;,  R)  =  .  mc‘,  i  =  1,2,... I! 

i  r  s  1 

i.e.  the  energy  deposited  to  the  i “  mode  by  each  electron 

traversing  the  cavity. 

c  :  initial  phase  angle  of  electron 

R  :  ring  on  which  electron  is  located  initially 
s 


DEFINITION  OF  EFFICIENCY  FOR  EACH  MODE 
,  N 


ni  = 


NW, 


r ,  s 


N:  Total  number  of  electrons 

W^:  Initial  energy  of  each  electron 

The  efficiency  for  each  mode  depends  on  the  amplitudes 
E01  Eq2,  ...»  EqM  of  all  the  modes,  i.e. 


ni  =  f i  (E01'  E02'  *  *  * '  EoM) 


i  =  1,  2,  ...»  M. 


TEST  OF  CONVERGENCE 

The  efficiency  for  each  mode  converges  to  a  fixed  value 
as  the  number  of  rings  increases.  (This  value  is  independent 
of  the  constant  phase  for  each  mode.) 
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Modes:  TE 


Oil'  TE012 
L  =  8  s  =  1 


r  =0.43 
o 

ii  =  3.5 

>  =1.14 

E.  =  0.1 

a  =1.5 

E2  =  0.05 

M-l 

c 

0) 

cn 

o  w 


u 

w 

^  a) 
c  £ 

•H 

U 

V£> 


Mode  TE 


>EST  OF  CONVERGENCE 


Mode  s : 

TE021 

TF 

1  “022 

L  =  5 

a  = 

1.5 

E1  =  0.1 

f 

?o  «  0.51 

s  = 

2 

E2  =  0.05 

.  =  1.14 

“1 

3.0 

♦ 

Number  of 

Rings 

Phase 

Difference 

n1(%) 

1 

0 

4.046 

1.071 

2 

0 

2.523 

0.587 

♦ 

3 

0 

2.512 

0.630 

4 

0 

2.513 

0.628 

4 

3.1415 

2.473 

0.630 

♦ 

5 

0 

2.493 

0.645 

5 

3.1415 

2.495 

0.644 

6 

0 

2.511 

0.629 

♦ 

6 

3.1415 

2.475 

0.631 

7 

0 

2.511 

0.629 

8 

0 

2.495 

0.625 

I 


DEFINITION  OF  STEADY  STATE 


Fi!E01/  E02 '  •*"  EoM  '  ?b}  "  r'iPb  Q^1 

1 

for  i  =  1 ,  2  ,  .  .  .  ,  M. 

r i?b  =  Pdi:  P°wer  deposited  to  the  i^'  mode  by  the  beam 


^  .  W  . 

— -Ar— — ■ —  =  P..  :  Wave  power  of  the  i^  mode 
y  •  j-i 


The  M  nonlinear  equations 

Fi(E0l’  E02’  Eo;.]'  Pb}  =  °'  1  ~  1'  2'  •**'  M 

define  the  steady  state  for  each  value  of  the  beam 

power  Pfa. 

At  the  steady  state,  there  is  a  balance  of  power 
for  each  mode. 
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BASIC  ASSUMPTION  OF  MODE  EXCITATION 

•  The  beam  voltage  is  fixed. 

•  As  the  beam  current  (or  Pfa)  is  increased  from 
zero  to  a  certain  value ,  the  system  traces  a 
sequence  of  steady  states  (as  defined  above) . 

•  The  mode  with  the  lowest  threshold  is  excited 
first . 

•  As  the  beam  current  is  increased,  a  new  mode  is 
excited  if  this  is  allowed  by  the  presence  of 
the  other  modes  in  the  cavity. 

•  The  above  assumption  is  valid  since: 

Rising  (or  decay)  time  of  fields  %  nsec 
Rising  time  of  beam  current  usee 

•  A  nonlinear  algebraic  system  may  have  more  than 

a  single  real  solution  for  given  Pfa,  and  only  one 
of  these  solutions  is  physically  acceptable. 
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Basic  results  will  be  presented  for  two  modes  as  plots  of  the  form 


•  Each  point  in  the  above  plots  is  obtained  by  solving  the  nonlinear 
algebraic  equations 


F1 (E01' 

Eo2'  V 

=  0 

F2  (E01' 

E02-  Pb> 

=  0 

for  given  P^. 


Iteration  algorithm  for  solving  the  M  nonlinear  algebraic  equations 


Fi(E01'  E02 '  EoK'  Pb}  •••<  M 


for  given  P^. 


Choose  arbitrary  Eq^,  E02'*’,EoK 


Compute  the  efficiencies  r,^. 


Compute  the  quantities  ,  F2,  F^ 


'^|Fil  <  G  ^\Yes  f - — — — v 

for  i=l,2,  . .  .  ,M  - ►{  End  computation] 


Update  the  amplitudes 
EoVEoi  +  5Eoi'  1  -  S- 


Go  to 


.  .  ,  M 


Example  2:  Two  modes  =  given  value 


S 2 "  Unstable  S.S 
S  ,  :  Stable  S.S. 


Contour  =  0  (F2  =  0)  divides  plane  into  two  regions  where  F^ 


(F2  >  0)  and  F±  <  0  (F2  <  0) 


•  When  F^  >  0  (F2  >  0)  the  mode  1  (2)  amplifies 
When  F^  <  0  (F2  <  0)  the  mode  1  (2)  decays. 


•  At  S^:  F2  <  0  -*■  Eq2  decays  to  S^ 

At  S2:  F|  =  0,  F2  =  0.  For  infinitesimal  E^  (noise)  F^  >  0 

amplifies  to  S^  (Unstable  steady  state) 

At  S3:  F-^  >  0  -*■  Eq1  amplifies  to  S^ 

At  S^:  F^  =  0,  F2  =  0.  For  infinitesimal  EQ2  (noise)  F2  <  0 

decays  back  to  S.  (Stable  steady  state) 


it  follows  that,  for 


From  the  definition  of  I'. .  , 
all  P^,  at  a  steady  state  £i.e. 


. . . ,  M)  the  relation 


Ki(E01'  E0 2 '  ’ ’ '  EoM) 


holds . 


0,  i 

1,  2 


1,2, 

.  .  ,  M 


Therefore,  for  a  given  cavity  (Q^  fixed)  the 
contours  i  =  1,  2,  ...,  M,  provide 

the  sequence  of  steady  states  in  the  cavity  as 
the  beam  power  increases,  if  that  is  also  the 
case  (i.e.  P^  increases)  as  we  move  along  these 
contours • 


Example  3:  Two  modes, 


K_  -  contour  plane 


contour  plane 


Assume : 


5 » 1 

D.5,  all  n 

bl  < 

Fb2  b3 

thr 

thr 

bl  < 

Pb2 

Easy  to  show  that  if  =  0,  then 
F2  =  n2PbQ2  1  '  Q^Q ” 

It  follows  (since  K2/(Q2/Q1)  >  1  for  all  K2)  that  for  infinitesimal 
Eq2,  the  relation  F2  <  0  holds.  Thus,  mode  2  will  never  be  excited 
as  the  beam  power  increases. 


Example  4 : 


0.6 


Oil 


The  contour  clots 


jroviae  m: orrr.attcr. 


on  the  threshold,  bean-,  cower  for  each  mode . 


In  the  linear  regime  ti.e.  E 


—  t  —  —  /  <-  f 


.  K  -  . 

i  fi  =  Dthr 


Ihus  P.  •  Q.  is  computed  from  the  value  of  X. 
£  1  1  1 


at  the  origin.  (Since  -  .  E~  .  ,  K . .  E‘.  in 

l  oi  fi,  oi 

the  linear  regime,  X.  . .  is  independent  of  E  .) 

i  /  ■•  in  ci 


The  relation 


p  “hr 

‘  bi 


l ,  ■-  in 

VQ1 


holds . 


•  For  two  modes: 


If  V^l  <  K2,  !in  “  ?bf  >  PSr 


If  VQ!  >  *2, tin  -  Pbf  "  Pbf 
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Modes : 


Q0/Q 


Modes:  TF 


Contour  |j  I  <  1 1  s 


Modes:  TE . . . ,  TE_ 
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Modes : 


Contour  Clot  of  efficiency  for  Mode 


Modes : 


Contour  Plot  of  f’ffic 


Modes 


Beam  Power  (kW) 


Mode  1 ,  Mode 


Beam  power  (kW) 


Mode  l ,  Mode 


Contour  Plot  of  X_(”W) 


Mode 


for  Modp 


Mode 


. 


o  o  o 

•  •  • 

o  ©  o 

m  <m  m 


(%)  AOU0TOTJJ3 


200  600  1000  1400  1800  2200 

Beam  power  (kW) 


Beam  power  (kW) 


(kW) 


VC,  =  70  keV,  \  /v  „  =  1.8,  L/r  =  12,  r  /r  =  0.48 
b  o  zo  w  o  w 

Constant  magnetic  field  profile 


operating  point 
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PRELIMINARY  CONCLUSIONS 
The  Q's  for  each  mode  are  important  in  the 
determination  of  mode  competition. 

For  modes  with  azimuthal  symmetry,  there  is 
indication  that,  as  the  beam  power  builds  up, 
the  mode  with  axial  eigennumber  <  =  1  will 
eventually  exclude  all  other  higher  order  modes 
in  the  cavity . 

There  is  indication  that  a  mode  with  infinite 
threshold  beam  power  may  be  excited  through  the 
excitation  of  other  modes  first  and  subsequent 


mode  conversion. 
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